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Abstract 

We study a large class of critical two-dimensional Ising models, namely critical 
Z -invariant Ising models. Fisher [Fis66| introduced a correspondence between 
the Ising model and the dimer model on a decorated graph, thus setting dimer 
techniques as a powerful tool for understanding the Ising model. In this paper, 
we give a full description of the dimer model corresponding to the critical Z- 
invariant Ising model, consisting of explicit expressions which only depend on the 
local geometry of the underlying isoradial graph. Our main result is an explicit local 
formula for the inverse Kasteleyn matrix, in the spirit of [Kcn02 , as a contour 
integral of the discrete exponential function of [Mer01a| IKen02j multiplied by a 
local function. Using results of [BdT08| and techniques of |dT07b| IKen02| , this 
yields an explicit local formula for a natural Gibbs measure, and a local formula 
for the free energy. As a corollary, we recover Baxter's formula for the free energy 
of the critical Z-invariant Ising model [Bax89| . and thus a new proof of it. The 
latter is equal, up to a constant, to the logarithm of the normalized determinant 
of the Laplacian obtained in [Ken02j . 



1 Introduction 

In |Fis66j . Fisher introduced a correspondence between the two-dimensional Ising model 
defined on a graph G, and the dimer model defined on a decorated version of this graph. 
Since then, dimer techniques have been a powerful tool for solving pertinent questions 
about the Ising model, see for example the paper of Kasteleyn |Kas67] . and the book 
of Mc Coy and Wu |MW73j . In this paper, we follow this approach to the Ising model. 

We consider a large class of critical Ising models, known as critical Z -invariant Ising 
models, introduced in |Bax86j . More precisely, we consider Ising models defined on 
graphs which have the property of having an isoradial embedding. We suppose that 

*Laboratoire de Probabilites et Modeles Aleatoires, Universite Paris VI Pierre et Marie Curie, 
Case courrier 188, 4 place Jussieu, F-75252 Paris CEDEX 05. cedric. boutillierSupmcfr. 
Supported in part by the Swiss National Foundation Grant 200020-120218/1. 

^Institut de Mathematiques, Universite de Neuchatel, Rue Emile-Argand 11, CII-2007 
Neuchatel. beatrice.detiliere@unine.ch. Supported in part by the Swiss National Founda- 
tions grants 47102009 and 200020-120218/1. 



the Ising coupling constants naturally depend on the geometry of the embedded graph, 
and are such that the model is invariant under star-triangle transformations of the 
underlying graph, i.e. such that the Ising model is Z -invariant. We suppose moreover 
that the coupling constants are critical by imposing a generalized self-duality property. 
The standard Ising model on the square, triangular and honeycomb lattice at the critical 
temperature are examples of critical Z-invariant Ising models. In the mathematics 
literature, this model has been studied by Mercat [MerOlbj . where the author proves 
equivalence between existence of Dirac spinors and criticality; by Chelkak and Smirnov 
[CS08j . where the authors prove conformal invariance in the scaling limit; in [BdTOS] . 
where we give a complete description of the equivalent dimer model, in the case where 
the underlying graph is periodic. The critical Z-invariant Ising model has also been 
widely studied in the physics literature, see for example |Bax86[ IAYP87[ IAP04t IAP071 
IMar97[[Ma?98l[CSn6] . 

Let G = {V{G), E{G)) be an infinite, locally finite, isoradial graph with critical coupling 
constants on the edges. Then by Fisher, the Ising model on G is in correspondence with 
the dimer model on a decorated graph with a well chosen positive weight function 
u on the edges. We refer to this model as the critical dimer model on the Fisher graph 
Q of G. It is defined as follows. A dimer configuration of ^ is a subset of edges of 
Q such that every vertex is incident to exactly one edge of M. Let M{Q) be the set 
of dimer configurations of Q. Dimer configurations of Q are chosen according to a 
probability measure, known as Gibbs measure, satisfying the following properties. If 
one fixes a perfect matching in an annular region of Q, then perfect matchings inside 
and outside of this annulus are independent. Moreover, the probability of occurrence 
of an interior matching is proportional to the product of its edges weights given by the 
weight function z/. 

The key objects used to obtain explicit expressions for relevant quantities of the dimer 
model are the Kasteleyn matrix, denoted by K, and its inverse. A Kasteleyn matrix 
is an oriented adjacency matrix of the graph Q, whose coefficients are weighted by the 
function v, introduced by Kasteleyn |Kas61j . Our main result is Theorem [H consisting 
of an explicit local expression for an inverse of the Kasteleyn matrix K. It can 
loosely be stated as follows, refer to Section [^2] for detailed definitions and to Theorem 
[S]for a precise statement. 

Theorem 1. Let x,y be two vertices of Q. Then the infinite matrix K^^ , whose 
coefficient K^y is given below, is an inverse Kasteleyn matrix. 

K~l = jf fx{X)fy{-X) Expx^y(A) log AdA + G^^y, 

where fx is a complex-valued function depending on the vertex x only; Exp,^ y is the 
discrete exponential function introduced in IMerOla!^ . see also IKenOS^ : Cx,y is a con- 
stant equal to when x and y are close, and else; Cx^y is a simple closed curve 
oriented counterclockwise containing all poles of the integrand, and avoiding a half-line 
dx,y starting from zero. 
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Before stating implications of Theorem [T] on the critical dimer model on the Fisher 
graph Q, let us make a few comments. 

1. Theorem [1] is in the spirit of the work of Kenyon |Ken02j , where the author obtains 
explicit local expressions for the critical Green's function defined on isoradial 
graphs, and for the inverse Kasteleyn matrix of the critical dimer model defined 
on bipartite isoradial graphs. Surprisingly both, the expression for the Green's 
function of |Ken02] and our expression for the inverse Kasteleyn matrix involve the 
discrete exponential function. This relation is pushed even further in Corollary 
\T2\ see the comment after the statement of Theorem [3l Note that the expression 
for the inverse Kasteleyn matrix obtained in |Ken02j does not hold in our setting, 
since the dimer model we consider is defined on the Fisher graph Q of G, which 
is not bipartite and not isoradial. 

2. Theorem [1] should also be compared with the explicit expression for the inverse 
Kasteleyn matrix obtained in |BdT08 ]. in the case where the graph Q is periodic. 
Then, by the uniqueness statement of Proposition 5 of [BdTOSj , the two expres- 
sions coincide. The proofs are nevertheless totally different in spirit, and we have 
not yet been able to understand the identity by an explicit computation, except 
in the case where G = J?. 

3. The most interesting features of Theorem [T] are the following. First, there is no 
periodicity assumption on the graph Q. Secondly, the expression for is local., 
meaning that it only depends on the local geometry of the underlying isoradial 
graph G. More precisely, it only depends on an edge-path of G between vertices 
X and y of G, naturally constructed from x and y. This implies that changing 
the isoradial graph G away from this path does not change the expression for 
K~^. Thirdly, explicit computations of become tractable, whereas, even in 
the periodic case, they remain very difficult with the explicit expression for 
given in |BdT08j . 

4. The structure of the proof of Theorem [1] is taken from jKen02] . The idea is to 
find complex-valued functions that are in the kernel of the Kasteleyn matrix K, 
then to define as a contour integral of these functions, and to define the 
contours of integration in such a way that KK~^ = Id. The great difficulty lies 
in actually finding the functions that are in the kernel of iiT, since there is no 
general method to construct them; and in defining the contours of integration. 
Indeed, the Fisher graph Q is obtained from an isoradial graph G, but has a more 
complicated structure, so that the geometric argument of |Ken02| does not work. 
The proof of Theorem [T] being long, it is postponed until the last section of this 
paper. 

Using Theorem [H an argument similar to [dT07b] , and the results obtained in |BdT08] . 
yields Theorem [2] giving an explicit expression for a Gibbs measure on Ai {Q) . When 
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the graph Q is periodic, this measure coincides with the Gibbs measure obtained in 
[BdT08] as weak hmit of Boltzmann measures on a natural toroidal exhaustion of the 
graph Q. A precise statement is given in Theorem [9] of Section [5.11 

Theorem 2. There is a unique Gibbs measure V defined on M{Q), such that the prob- 
ability of occurrence of a subset of edges {ei = xiyi, • • • , XkVk} in a dimer configuration 
of Q, chosen with respect to the Gibbs measure V is: 

where is given by Theorem{^ and {K~^)[xi,yi,-- ,Xk,yk} ^"^ sub-matrix of , 

whose lines and columns are indexed by vertices • • • ,Xk,yk}- 



Theorem [2] is a result about Gibbs measures with no periodicity assumption on the 
underlying graph. There are only very few examples of such instances in statistical me- 
chanics. Moreover, the Gibbs measure V inherits the properties of the inverse Kasteleyn 
matrix of Theorem (T) it is local, and allows for explicit computations, examples 
of which are given in Appendix 1X1 

Let us now assume that the graph G is Z^-periodic. Using Theorem [H and the tech- 
niques of |Ken02| . we obtain an explicit expression for the free energy of the critical 
dimer model on the graph Q, see Theorem 1101 depending only on the angles of the fun- 
damental domain. Using Fisher's correspondence between the Ising and dimer models, 
this yields a new proof of Baxter's formula for the free energy of the critical Z-invariant 
Ising model, denoted by //, see Section [52] for definitions. 

Theorem 3 f [Bax89j ). 

/, = -|y(Gi)|i^- V [^logtan0e + - toe)+i(5-0. 

e&E{Gi) 



Note that the free energy // of the critical Z- invariant Ising model is, up to a multiplica- 
tive constant — ^ and an additive constant, the logarithm of the normalized determinant 
of the Laplacian obtained by Kenyon |Ken02j . See also Corollary I12[ where we explic- 
itly determine the constant of proportionality relating the characteristic polynomials 
of the critical Laplacian and of the critical dimer model on the Fisher graph Q, whose 
existence had been estabhshed in [BdTOSj . 

Outline of the paper 



Section 2: Definition of the critical Z-invariant Ising model. 

Section 3: Fisher's correspondence between the Ising and dimer models. 
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Section 4: Definition of tlie Kasteleyn matrix K. Statement of Theorem [5l giving an ex- 
plicit local expression for the coefficient K~y of the inverse Kasteleyn matrix. 
Asymptotic expansion of K~y, as |a; — ?/| oo, using techniques of |Ken02j . 

Section 5: Implications of Theorem [S] on the critical dimer model on the Fisher graph Q: 
Theorem [9] gives an explicit local expression for a natural Gibbs measure, and 
Theorem [10] gives an explicit local expression for the free energy. Corollary: 
Baxter's formula for the free energy of the critical Z-invariant Ising model. 

Section 6: Proof of Theorem [5] 

Acknowledgments: we would like to thank Richard Kenyon for asking the questions 
solved in this paper. 



2 The critical Z-invariant Ising model 

Consider an unoriented finite graph G = {V{G), E{G)), together with a collection of 
positive real numbers J = {Je)e^E(G) indexed by the edges of G. The Ising model on 
G with coupling constants J is defined as follows. A spin configuration o" of G is a 
function of the vertices of G with values in {—1, +1}. The probability of occurrence of 
a spin configuration a is given by the Ising Boltzmann measure^ denoted P^: 

P\o) = exp I ^ JeCTufyv I , 

\e=uv(^E{G) / 

where Z'^ = J2ae{-i ij^C?) exp ^X^e=™e-B(G) Jef^uO'vj , is the Ising partition function. 

We consider Ising models defined on a class of embedded graphs which have an addi- 
tional property called isoradiality. A graph G is said to be isoradial [Ken02j . if it has an 
embedding in the plane such that every face is inscribed in a circle of radius 1, and all 
circumcenters of the faces are in the closure of the faces. Prom now on, when we speak 
of the graph G, we mean the graph together with a particular isoradial embedding in 
the plane. Examples of isoradial graphs are the square and the honeycomb lattice, see 
Pigure [T] (left) for a more general example of isoradial graph. 

To such a graph is naturally associated the diamond graph, denoted by G*, defined as 
follows. Vertices of G^ consist in the vertices of G, and the circumcenters of the faces 
of G. The circumcenter of each face is then joined to all vertices which are on the 
boundary of this face, see Pigure [T] (center). Since G is isoradial, all faces of are 
side-length-1 rhombi. Moreover, each edge e of G is the diagonal of exactly one rhombus 
of G*; we let Be be the half- angle of the rhombus at the vertex it has in common with 
e, see Pigure [D (right). 
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Figure 1: Left: example of isoradial graph. Center: corresponding diamond graph. 
Right: rhombus half-angle associated to an edge e of the graph. 



The same construction can be done for infinite and toroidal isoradial graphs, in which 
case the embedding is on a torus. When the isoradial graph is infinite and non periodic, 
in order to ensure that the embedding is locally finite, we assume that there exists an 
e > such that the half-angle of every rhombus of lies between e and § — This 
implies in particular that vertices of G have bounded degree. 

It is then natural to choose the coupling constants J of the Ising model defined on an 
isoradial graph G, to depend on the geometry of the embedded graph: let us assume 
that Je is a function of Oe-, the rhombus half- angle assigned to the edge e. 

We impose two more conditions on the coupling constants. First, we ask that the 
Ising model on G with coupling constants J as above is Z -invariant, that is, invariant 
under star-triangle transformations of the underlying graph. Next, we impose that the 
Ising model satisfies a generalized form of self- duality. These conditions completely 
determine the coupling constants J, known as critical coupling constants: for every 
edge e of G, 



The Z-invariant Ising model on an isoradial graph with this particular choice of cou- 
pling constants is referred to as the critical Z -invariant Ising model. This model was 
introduced by Baxter in |Bax86] . A more detailed definition is given in [BdT08] . 

3 Fisher's correspondence between the Ising and dimer 
models 

Fisher |Fis66j exhibits a correspondence between the Ising model on any graph G drawn 
on a surface without boundary, and the dimer model on a "decorated" version of G. 
Before explaining this correspondence, let us first recall the definition of the dimer 
model. 




(1) 
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3.1 Dimer model 



Consider a finite graph Q = {V{Q), E{Q)), and suppose that edges of Q are assigned a 
positive weight function u = {ve)e&E[g)- The dimer model on Q with weight function v 
is defined as follows. 

A dimer configuration MofQ, also called perfect matching, is a subset of edges of Q such 
that every vertex is incident to exactly one edge of M. Let Ai{G) be the set of dimer 
configurations of the graph Q. The probability of occurrence of a dimer configuration 
M is given by the dimer Boltzmann measure, denoted V^: 



where Z'^ = J^MeMiG) TleeM dimer partition function. 



3.2 Fisher's correspondence 

Consider an Ising model on a finite graph G embedded on a surface without boundary, 
with coupling constants J. We use the following slight variation of Fisher's correspon- 
dence |Fis66] ■ 

The decorated graph, on which the dimer configurations live, is constructed from G 
as follows. Every vertex of degree A; of G is replaced by a decoration consisting of 3k 
vertices: a triangle is attached to every edge incident to this vertex, and these triangles 
are linked by edges in a circular way, see Figure [2] below. This new graph, denoted by 
Q, is also embedded on the surface without boundary and has vertices of degree 3. It 
is referred to as the Fisher graph of G. 




Figure 2: Left: a vertex of G with its incoming edges. Right: corresponding decoration 
in Q. 

Fisher's correspondence uses the high temperature expansior0 of the Ising partition 

'^It is also possible to use the low temperature expansion, if the spins of the the Ising model do not 
sit on the vertices but on the faces of G. 
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function, see for example [Bax89j : 



Z'=[ H cosh(Je) 2l^(^)l Y^YltanhiJ,), 

\e(^E{G) / CePeeC 

where P is the family of all polygonal contours drawn on G, for which every edge of 
G is used at most once. This expansion defines a measure on the set of polygonal 
contours P of G: the probability of occurrence of a polygonal contour C is proportional 
to the product of the weights of the edges it contains, where the weight of an edge e is 
tanh( Je)- 

Here comes the correspondence: to any contour configuration C coming from the high- 
temperature expansion of the Ising model on G, we associate 2l^*^'^)l dimer configu- 
rations on Q: edges present (resp. absent) in C are absent (resp. present) in the 
corresponding dimer configuration of Q. Once the state of these edges is fixed, there 
is, for every decorated vertex, exactly two ways to complete the configuration into a 
dimer configuration. Figure [3] below gives an example in the case where G is the square 
lattice Z2. 



1 2 



3 4 




Figure 3: Polygonal contour of Z^, and corresponding dimer configurations of the 
associated Fisher graph. 



Let us assign, to an edge e oi Q, weight = 1, if it belongs to a decoration; and weight 
= coth Jg, if it corresponds to an edge of G. Then the correspondence is measure- 
preserving: every contour configuration C has the same number [2^^^^^^) of images by 
this correspondence, and the product of the weights of the edges in C, J^ggj- tanh( Jg) 
is proportional to the weight He^C coth(Je) of any of its corresponding dimer configu- 
rations for a proportionality factor, HeeECG) taiih(Je), which is independent of C. 

As a consequence of Fisher's correspondence, we have the following relation between 
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the Ising and dimer partition functions: 




H sinh(Je) 



(2) 



Fisher's correspondence between Ising contour configurations and dimer configurations 
naturahy extends to the case where G is an infinite planar graph. 

3.3 Critical dimer model on Fisher graphs 

Consider a critical Z-invariant Ising model on an isoradial graph G on the torus, or 
on the whole plane. Then, the dimer weights of the corresponding dimer model on the 
Fisher graph Q are: 



We refer to these weights as critical dimer weights, and to the corresponding dimer 
model as critical dimer model on the Fisher graph Q. 

4 Kasteleyn matrix on critical infinite Fisher graphs 

In the whole of this section, we let G be an infinite isoradial graph, and Q be the 
corresponding Fisher graph. We suppose that edges of Q are assigned the dimer critical 
weight function denoted by v. Recall that G^ denotes the diamond graph associated 
to G. 

4.1 Kasteleyn and inverse Kasteleyn matrix 

The key object used to obtain explicit expressions for the dimer model on the Fisher 
graph Q is the Kasteleyn matrix introduced by Kasteleyn in |Kas61j . It is a weighted, 
oriented adjacency matrix of the graph Q defined as follows. 

A Kasteleyn orientation of Q is an orientation of the edges of Q such that all elementary 
cycles are clockwise odd, i.e. when traveling clockwise around the edges of any elemen- 
tary cycle of Q, the number of co-oriented edges is odd. When the graph is planar, such 
an orientation always exists |Kas67| . For later purposes, we need to keep track of the 
orientation of the edges of Q. We thus choose a specific Kasteleyn orientation of Q in 
which every triangle of every decoration is oriented clockwise. Having a Kasteleyn ori- 
entation of the graph Q then amounts to finding a Kasteleyn orientation of the planar 
graph obtained from Q by contracting each triangle to a single vertex, which exists by 




1 if e belongs to a decoration, 

u{6e) = cot (^) if e comes from an edge of G. 
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Kasteleyn's theorem [Kas67) . Refer to Figure S] for an example of such an orientation 
in the case where G = Z^. 




Figure 4: An example of Kasteleyn orientation of the Fisher graph of Z^, in which 
every triangle of every decoration is oriented clockwise. 



The Kasteleyn matrix corresponding to such an orientation is an infinite matrix, whose 
rows and columns are indexed by vertices of Q, defined by: 

where 



1 if X ~ y, and x — > y 
— 1 if X ~ y, and x <— y 
else. 



Note that K can be interpreted as an operator acting on C^^^^: 

yev{g) 

An inverse of the Kasteleyn matrix K, denoted K~~^, is an infinite matrix whose rows 
and columns are indexed by vertices of Q, and which satisfies KK~^ = Id. 

4.2 Local formula for an inverse Kasteleyn matrix 

In this section, we state Theorem [S] proving an explicit local expression for the coef- 
ficients of an inverse K~^ of the Kasteleyn matrix K. This inverse is the key object 
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for the critical dimer model on the Fisher graph Q. Indeed, it yields an explicit local 
expression for a Gibbs measure on dimer configurations of Q, see Section [5.11 It also 
allows for a simple derivation of Baxter's formula for the free energy of the critical 
Z-invariant Ising model, see Section [5.21 

This section is organized as follows. The explicit expression for the coefficients of 
given by Theorem below, see also Theorem [H is a contour integral of an integrand 
involving two quantities: 

• a complex- valued function depending on a complex parameter, and on the vertices 
of the graph Q only, defined in Section 14.2.21 

• the discrete exponential function which first appeared in [MerOlaj , see also |Ken02j . 
It is a complex valued function depending on a complex parameter, and on an 
edge-path between pairs of vertices of the graph G, defined in Section I4.2.3[ 

Theorem [S] is then stated in Section [4.2.41 Since the proof is long, it is postponed until 
Section m In Section [4.2.5l using the same technique as |Ken02| . we give the asymptotic 
expansion of the coefficient K~y of the inverse Kasteleyn matrix, as |x — y| — > oo. 

4.2.1 Preliminary notations 

From now on, vertices of Q are written in normal symbol, and vertices of G in boldface. 
Let X be a vertex of G, then x belongs to the decoration corresponding to a unique 
vertex of G, denoted by x. Conversely, vertices of G of the decoration corresponding 
to a vertex x of G are labeled as follows, refer to Figure [5] for an example. Let d{x) 
be the degree of the vertex x in G, then the corresponding decoration of G consists of 
d{x) triangles, labeled from 1 to d{x) in counterclockwise order. For the /c-th triangle, 
let Vkix) be the vertex incident to an edge of G, and let Wfc(x), Zfc(x) be the two other 
vertices, in counterclockwise order, starting from Vfe(x). 



Later on, when no confusion occurs, we will drop the argument x in the above labeling. 
Define a vertex a; of ^ to be of type 'v if a; = Vfc(x) for some k £ {1, • • • , d{x)}, and 
similarly for 'u)' and 




G 




Figure 5: Notations for vertices of G- 
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The isoradial embedding of the graph G fixes an embedding of the corresponding di- 
amond graph . There is a natural way of assigning rhombus unit- vectors of to 
vertices of Q: for every vertex x of G, and every k G {!,•■■ let us associate 

the rhombus unit- vector e°'^k^^) to Wk{^), 6*°^*='''' to Zk{x), and the two rhombus- unit 
vectors e^^'^fcW, e*°^fe('') to Vfc(x), as in Figure El below. Note that e*°"'fe('') = e^^^+ii^) , 




Figure 6: Rhombus vectors of the diamond graph G^ assigned to vertices of Q. 



4.2.2 Complex-valued function on the vertices of Q 

Let us introduce the complex-valued function defined on vertices of Q and depending 
on a complex parameter, involved in the integrand of the contour integral of given 
by Theorem O Define / : V{g) x C ^ C, by: 



/(wfc(x), A) := /^fe(x)(A) 

/(zfc(x), A) := /2fc(x)(A) 

/(wfc(x), A) := /t,fc(x)(A) 

for every x £ G, and every k € {1, • • • ,d(x)}. In order for the function / to be well 
defined, the angles awf,(x)j '^z^ix) need to be well defined mod 47r, indeed half-angles 
need to be well defined mod 2tt. Let us define them inductively as follows, see also 
Figure El The definition strongly depends on the Kasteleyn orientation introducted in 
Section ICT Fix a vertex xq of G, and set a^i(xo) = 0- Then, for vertices of G in the 
decoration of a vertex x G G, define: 

Q^Zfc(x) + 2^A;(x), where Ok{x) > is the rhombus half- angle of Figure [71 

{dwkix) if ttie edge it;fc(x)zfc+i(x) is oriented from Wkix) to Zfc+i(x) 

««>fc(x) + 2vr else. 

(3) 



e 2 



e 2 



= /«lfe(x)(A) + /2fe(x)(A), 
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Here is the rule defining angles in the neighboring decoration, corresponding to a vertex 
y of G. Let k and £ be indices such that ffe(x) is adjacent to vi{y) in Q. Then, define: 




Figure 7: Notations for the definition of the angles in R/47rZ. 

Lemma 4. For every vertex x of G, and every k G {1, • • • ,d(x)}, the angles a^j.(x)! 
'^2ft(x)) well defined in ]R/47rZ. 

Proof. It suffices to show that when doing the inductive procedure around a cycle of G, 
we obtain the same angle modulo 47r. There are two types of cycles to consider: inner 
cycles of decorations, and cycles of Q coming from the boundary of a face of G. These 
cycles are represented in Figure [HI 




Figure 8: The two types of cycles on which the definition of the angles in M/47rZ needs 
to be checked. Left: Ci, the inner cycle of a decoration, oriented counterclockwise. 
Right: C2, a cycle coming from the boundary of a face of G, oriented clockwise. 



Let Ci be a cycle of the first type, that is Ci is the inner cycle of a decoration corre- 
sponding to a vertex x of G, oriented counterclockwise. Let d = c?(x) be the degree of 
the vertex x in G, then: 

Ci = {zi{x),wi{x), . . . , Zd{x),Wdi:x.), zi{x)) = {zi,wi,...,Zd,Wd,zi). 



13 



According to our choice of Kasteleyn orientation, all edges of this cycle are oriented 
counterclockwise with respect to the face surrounded by Ci, except an odd number n 
of edges of the form w^Zk+i- By definition of the angles ([3|), when jumping from to 
Zk+i, the change in angle is 20^ if WkZk+i is oriented counterclockwise, and 26k + 27r 
otherwise. The total change along the cycle is thus, 

d 

2ek + 27rn = 27r(n + 1) = mod An. 

k=l 

Let C2 be a cycle of the second type, around a face touching m decorations correspond- 
ing to vertices xi, • • • ,Xm of G. Suppose that C2 is oriented clockwise. 

C2 = (u;^i(xi),2;£^+i(xi),i;£^+i(xi ),..., 

• • ■ ,Ve^{Xm),We^ {Xm),Ze^+i{Xm),Ve^+i{Xm),Ve^{xi),Wi^ (xi)). 

The cycle C2 contains 4m edges, and the number n of co-oriented edges along C2 is odd 
by definition of a Kasteleyn orientation. By our choice of Kasteleyn orientation, the only 
edges that may be co-oriented along the cycle, are either of the form w^^. (x,fc)2;£^_(_i(xfc), 
or of the form Vi^+i{xk)vi^_^^{xk+i). 

For every k G {!,••• there are two different contributions to the total change in 
angle from w^^{xk) to u;^^-^ (xfc+i). First, from we^{xk) to ze^+i{xk), by there is 
a contribution of 27r if the corresponding edge is co-oriented. Then, from 2:^^_|_i(xfc) to 
w;4^-^(xfc+i), by (jl]), the angle is changed by 29i^^+i{xk) - n = 26*4^^ (xfc+i) - vr, plus an 
extra contribution of 27r if the edge V£^,^i{x)vii^^^{xk+i) is co-oriented. Thus the total 
change in angle is: 

m 

J];(204(xfe) -7r) + 27rn. 
k=l 

But the terms vr — 29i^{xk) are the angles of the rhombi at the center of the face 
surrounded by C2, which sum to 27r. Therefore, the total change is equal to 27r(n — 1) 
which is congruent to mod An. □ 

These angles at vertices are related to the notion of spin structure on surface graphs 
(see |CR07j . [CR08j . |Kup98| ). A spin structure on a surface is equivalent to the data 
of a vector field with even index singularities. Kuperberg explained in |Kup98| how to 
construct from a Kasteleyn orientation a vector field with odd index singularities at 
vertices of G. One can then obtain the spin structure by merging the singularities into 
pairs using a reference dimer configuration. 

The angles a«,j,(x) and ^^^.(x) defined here are directly related to the direction of the vec- 
tor field with even index singularities obtained from Kuperberg's construction applied 
to our choice of Kasteleyn orientation, and from the pairing of singularities correspond- 
ing to the following reference dimer configuration: 
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• ffc(x) <-> f^(y) if they are neighbors. 
4.2.3 Discrete exponential functions 

Let us define the discrete exponential function, denoted Exp, involved in the inte- 
grand of the contour integral of given by Theorem [5l This function first ap- 
peared in [MerOlaj . see also [Ken02| . In order to simplify notations, we use a dif- 
ferent labeling of the rhombus vectors of G^. Let x, y be two vertices of G, and let 
y = xi,X2,--- ,x„4.i = X be an edge-path of G from y to x. The complex vector 
Xj+i — Xj is the sum of two unit complex numbers e*^^ + e^'^^ representing edges of the 
rhombus in associated to the edge XjXj+i. 

Then, Exp : V{G) x V{G) x C ^ C is defined by: 

Exp(.,y, A) := Exp.,(A) = (sj±|) {^"j . 

j = l 

The function is well defined (independent of the choice of edge-path of G from y to x) 
since the product of the multipliers around a rhombus is 1. 



4.2.4 Inverse Kasteleyn matrix 



We now state Theorem [5] proving an explicit local formula for the coefficients of an 
inverse of the Kasteleyn matrix K. The vertices x and y of ^ in the statement 
should be thought of as being one of tffc(x), Zfc(x), ffc(x) for some x G G and some 
k G {!,••• ,(i(x)}, and similarly for y. The proof of Theorem [5] is postponed until 
Section [H 

Theorem 5. Let x,y be any two vertices of Q . Then the infinite matrix K^^, whose 
coefficient K~y is given by ([5]) below, is an inverse Kasteleyn matrix. 

^ ^ /.(A)/^(-A)Exp,,y(A)logAdA + a,j,. (5) 



The contour of integration Cx,y is a simple closed curve oriented counterclockwise con- 
taining all poles of the integrand, and avoiding the half-line d^^y starting from zero H. 
The constant Gx,y is given by 



Gx' 



1 

4 

_i 

4 

(-1) 



n(x.y) 







ifx = 
ifx = 
i/x = 
else, 



y = Wk{x) 
y = ^fc(x) 

y, x ^ y, and x and y are of type 'w 'or 'z' 



^ In most cases, the half-line dx,y is oriented from x to y, where x, y are vertices of at distance at 
most two from x, y, constructed in Section l673l Refer to this section and Remark [15] for more details. 
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where n{x, y) is the number of edges oriented clockwise in the clockwise arc from x to 
y of the inner cycle of the decoration corresponding to x. 

Remark 6. The explicit expression for K~y given in ^ has the very interesting feature 
of being local, i.e. it only depends on the geometry of the embedding of the isoradial 
graph G on a path between x and y, where x and y are vertices of at distance at 
most 2 from x and y, constructed in Section 16.31 A nice consequence of this property 
is the following. For i = 1,2, let Gi be an isoradial graph with corresponding Fisher 
graph Qi. Let Ki be the Kasteleyn matrix of the graph Qi, whose edges are assigned 
the dimer critical weight function. Let be the inverse of Ki given by Theorem 

O If G\ and G2 coincide on a ball B, and if the Kasteleyn orientations on Qi and Q2 
are chosen to be the same in the ball B, then for every couple of vertices (x, y) in B at 
distance at least 2 from the boundary, we have 



4.2.5 Asymptotic expansion of the inverse Kasteleyn matrix 



As a corollary to Theorem[5l and using computations analogous to those of Theorem 4.3 
of |Ken02j . we obtain the asymptotic expansion for the coefficient K~y of the inverse 
Kasteleyn matrix, as |x — y| — > 00. In order to give a concise statement, let us introduce 
the following simplified notations. When x is of type 'w' or 'z', let e*" denote the 
corresponding rhombus unit-vector of G*; and when x is of type let e*"^, e*"^ be 
the two corresponding rhombus unit-vectors of G*, defined in Section [4. 2. 1[ Moreover, 
set: 

{1 if j; is of type 'w' or 'f ' 
— 1 if a; is of type 'z'. 

A superscript "prime" is added to these notations for the vertex y, and ey is defined in 
a similar way. 

Corollary 7. The asymptotic expansion, as |x — y| — > 00, of the coefficient K~y of the 
inverse Kasteleyn matrix of Theorem\^ is: 



K 



x,y 



27r 



x-y 



+ 



|x~y| 



|x-y| 



if X and y are of type 'w ' or 'z 

if X is of type 'w ' or 'z \ y is of type 'v ', 

if X and y are of type 'v '. 



Proof. We follow the proof of Theorem 4.3 of |Ken02j . By Theorem O y is given by 
the integral: 



K-l = [ f.{X)fy{-X) Exp,,y(A) log AdA + 



(6) 
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where Cx^y is a simple closed curve oriented counterclockwise containing all poles of the 
integrand, and avoiding the half-line dx,y starting from in the direction from x to 
y, where x, y are vertices of G^, at distance at most 2 from x and y, constructed in 
Section [6^31 We choose dx,y to be the origin for the angles, i.e. the positive real axis 
R+. The coefficient Cx,y equals to whenever x and y are not in the same decoration. 
It will be a fortiori the case when they are far away from each other. 

It suffices to handle the case where both x and y are of type 'w' or 'z'. The other cases 
are then derived using the relation between the functions /: 

/...(A) = /«,,(A) + /.,(A). 

Define, 

n 

^ = x-y = ^e^^^ +e^^^-. 

where e*^^ , e*"^^ , • • • , e*^" , e*'''" are the steps of a path in from y to x. When x and 
y are far apart, since x and y are at distance at most 2 from x and y, the direction of 
xy is not very different from that of xy, ensuring that 3fi(^) < 0. 

The contour Cx^y can be deformed to a curve running counterclockwise around the ball 
of radius R {R large) around the origin from the angle to 2tt, then along the positive 
real axis, from R to r {r small), then clockwise around the ball of radius r from the 
angle 2tt to 0, and then back along the real axis from r to R. The rational fraction 
/2:(A)/j,(— A) Expx y(A) behaves like 0(1) when A is small, and like O (p-), when A is 
large. As a consequence, the contribution to the integral around the balls of radius r 
and R converges to 0, as we let r — > and — > oo. The logarithm differs by 2i7r on 
the two sides of the ray M+, the integral ([6]) is therefore equal to: 

1 roo 



As in |Ken02j . when |x — y| is large, the main contribution to this integral comes from 
a neighborhood of the origin and a neighborhood of infinity. When A is small, we have: 

e-t+o(A)^ 4l- = e-^i^ow^ ^ = expf2e-/5A + 0(A3) 



e*" - A ' e*" + A ' e^>^ - X 
Thus, for small values of A: 

e*f -A- e^^J + A e*'^^' + A 

fx{X)fy{-X) Exp^JX) = cxey ^.^ _ ^ ^^ [[ ^.^^ _ ^ ^^^^ _ ^ 

= €xeye-'^ exp (2^A + 0(A) + O(nA^)) . 
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Integrating this estimate from to gives 

1 1 
J^^ f,Wfyi-X)Exp^^y{X)dX = €,ey j^^ e-'^exp (2C"A + 0(n-^)) dA 



Similarly, when A is large, we have: 

e^" - A " A ' e^°' + A " A ' e^^ - A 

Thus for large values of A, 



exp (2e^^A-^ + 0(A- 



/,(A)/j,(-A)Exp,,y(A) = -e,ej^^^exp {2^\-^ + 0{\-^) + 0{n\-^)) . 
Computing the integral of this estimate for A G {\/n, oo) gives 



j-OO 

I fx{X)fy{-\) Expx_y(A)dA = e^ey 



The rest of the integral between n~2 and na is negligible (see |Ken02j ) . As a conse- 
quence, 

□ 



Remark 8. Using the same method with an expansion of the integrand to a higher 



order would lead to a more precise asymptotic expansion of K^l. 



5 Critical dimer model on infinite Fisher graphs 

Let Q be an infinite Fisher graph obtained from an infinite isoradial graph G. Assume 
that edges of Q are assigned the dimer critical weight function u. In this section, we give 
a full description of the critical dimer model on the Fisher graph Q, consisting of explicit 
expressions which only depend on the local geometry of the underlying isoradial graph 
G. More precisely, in Section 15.11 using the method of |dTn7bj . we give an explicit 
local formula for a Gibbs measure on dimer configurations of Q, involving the inverse 
Kasteleyn matrix of Theorem [3 When the graph is periodic, this measure coincides 
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with the Gibbs measure of |BdT08| , obtained as weak hmit of Boltzmann measures on 
a natural toroidal exhaustion. Then, in Section 15.21 we assume that the graph Q is 
periodic and, using the method of |Ken02j . we give an explicit local formula for the free 
energy of the critical dimer model on Q. As a corollary, we obtain Baxter's celebrated 
formula for the free energy of the critical Z-invariant Ising model. 



5.1 Local formula for the critical dimer Gibbs measure 

A Gibbs measure on the set of dimer configurations of G, is a probability measure 

on M{Q), which satisfies the following. If one fixes a perfect matching in an annular 
region of Q, then perfect matchings inside and outside of this annulus are independent. 
Moreover, the probability of occurrence of an interior matching is proportional to the 
product of the edge weights. 

In order to state Theorem [9l we need the following definition. Let if^ be a finite subset 
of edges of Q. The cylinder set is defined to be the set of dimer configurations of 
Q containing the subset of edges £; it is often convenient to identify £ with Ag. Let 
be the cr-algebra of A4{G) generated by the cylinder sets {As)^^^^-^^- Recall that 
K denotes the infinite Kasteleyn matrix of the graph Q, and let be the matrix 
inverse given by Theorem El 

Theorem 9. There is a unique probability measure V on {M{Q),T), such that for 
every finite collection of edges £ = {ei = xiyi, • • • , = x^yk} C E{Q), the probability 
of the corresponding cylinder set is 

V{Ae)=V{e^,--- ,ek) = (U^-^^v^ Pf ( , (7) 

where K^^ is given by Theorem\^ and ,a;fc,j/fc} sub-matrix of K^^ 

whose rows and columns are indexed by vertices {xi,yi, ■ ■ ■ ,Xk,yk}- 

Moreover, V is a Gibbs measure. When Q is 1? -periodic, V is the Gibbs measure 
obtained as weak limit of the Boltzmann measures Vn on the toroidal exhaustion {Qn = 
G/{nI?)]n>i ofG. 



Proof. The idea is to use Kolmogorov's extension theorem. The structure of the proof 
is taken from |dT07bj . Let iS = {ei = xiyi, . . . = Xnyn} be a finite subset of edges 
of Q. Denote by the cr-algebra generated by the cylinders [Ag/)^,^^. We define 
an additive function Vs on J^^, by giving its value on every cylinder set A£r, with 

£ {^il 5 • • • ! Sjfc } C iS. 
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It is a priori not obvious that Vs defines a probability measure on T£ . Let us show that 
this is indeed the case. Let be the subset of vertices of G to which the decorations 
containing xi,yi, . . . , Xn, Un retract. Define Q to be a simply-connected subset of rhombi 
of containing all rhombi adjacent to vertices of V^. 

By Proposition 1 of |dT07b] . there exists a Z^-periodic rhombus tiling of the plane con- 
taining Q, which we denote by G^^. Moreover, the rhombus tiling G*^ is the diamond 
graph of a unique isoradial graph G^ whose vertices contains the subset V^. Let 
be the Fisher graph of G^, then the graphs G and coincide on a ball B containing 

2^1 ) yi ) • • • ) ^ni Un- 

Endow edges of with the critical weights and a periodic Kasteleyn orientation, 
coinciding with the orientation of the edges of Q on the common ball B. We know by 
Remark El that the coefficients of the inverses {KP)~^ and given by Theorem 
are equal for all pairs of vertices in B. 

On the other hand, we know by Corollary [7| that the coefficient {KP)~^ 0, as 
|x — y| — > cxD. By Proposition 5 of |BdT08| . stating uniqueness of the inverse Kasteleyn 
matrix decreasing at infinity in the periodic case, we deduce that {KP)~^ is in fact the 
inverse computed in [ BdTOSj by Fourier transform. In Theorem 6 of |BdT08j . we use 
the inverse (KP)~^ to construct the Gibbs measure on obtained as weak 

limit of Boltzmann measures on the natural toroidal exhaustion of Q^, which has the 
following explicit expression: let {e'l = x'lu'i, - ■ ■ , e'^ = x'^y^} be a subset of edges of Q^, 
then 



In particular, the expression of the restriction of to events involving edges in £ 
(seen as edges of Q^) is equal to the formula defining Vs- As a consequence, Vs is a 
probability measure on Tg. 

The fact that Kolmogorov's consistency relations are satisfied by the collection of prob- 
ability measures (Vs)^ is immediate once we notice that any finite number of these 
probability measures can be interpreted as finite-dimensional marginals of a Gibbs 
measure on dimer configurations of a periodic graph with a large enough fundamental 
domain. The measure V of Theorem [9] is thus the one given by Kolmogorov's extension 
theorem applied to the collection (Vs)^. The Gibbs property follows from the Gibbs 
property of the probability measures V^. □ 

The probability of single edges are computed in details in Appendixj^ using the explicit 
form for given by Theorem O Consider a rhombus of with vertices x,t,y,u: 
x and y are vertices of G; t and u are vertices of G* . Let 6 be the half-angle of this 
rhombus, measured at x. In Q, we have that ii v = Vfc(x) and v' = vi^y) are end points 
of the edge e coming from xy in the decoration process, then we have 



^^(e;,...eD= ^<;JPfC(^V;,....a)■ 




71-26 
271 cos 9 



(8) 
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In Fisher's correspondence, the presence of the edge e in the dimer configuration corre- 
sponds to the fact that xy is not covered by a piece of contour. In the low temperature 
expansion of the Ising model, this thus means that the two spins at t and u have the 
same sign. Since they can be both + or — , the probability that they are both +, as a 
function of (f) = ^ — 9, the angle of the rhombus measured at u, is: 

1/1 vr - 26* \ _ 1 
2 V2 27r cos 9 7 4 27rsin>' 




5.2 Free energy of the critical dimer model 

In this section, we suppose that the isoradial graph G is periodic, so that the corre- 
sponding Fisher graph Q is. A natural exhaustion of Q by toroidal graphs is given by 
{Gn}n>ij where Qn = Q/n7?, and similarly for G. The graphs Qi and Gi are known 
as the fundamental domains of G and G respectively. In order to shorten notations in 
Theorem [To] below and in the proof, we write \E„\ = and = 

The free energy per fundamental domain of the critical dimer model on the Fisher graph 
G is denoted /d, and is defined by: 

fD = - hm ^logz;;, 

n—>oo n 

where is the dimer partition function of the graph Gn, whose edges are assigned the 
critical weight function i^. 

Theorem 10. The free energy per fundamental domain of the critical dimer model on 
the Fisher graph G is given by: 

/d = -(|£iI+|1'i|)!^ + 



^-26*6, , 1, 9e 1 ^ , /vr 

— — — log tan Oe- - log cot y - - \^L{9e) + L 

where L is the Lobachev sky function, L{x) = — Jq^ log 2 sin i dL 



Proof. For the proof, we follow the argument given by Kenyon [Ken02j to compute the 
normalized determinant of the Laplacian and the Dirac operator on isoradial graphs. 
See also |dT07aj for a detailed computation of the free energy of dimer models on 
isoradial bipartite graphs. 

In |BdT08j . we proved that 

„ Iff. , />/ . dz dw 

JD = -- logdet K{z,w) 



2 7 Jj2 2i'Kz 2i'Kw 
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where K{z, w) is the Fourier transform of the infinite periodic Kasteleyn matrix of the 
critical dimer model on Q. 



We need the following definition of |KS05| . Recall that is the diamond graph 
associated to the isoradial graph G. A train-track of is a path of edge-adjacent 
rhombi of G^, which does not tmn: on entering a face, it exits along the opposite edge. 
As a consequence, each rhombus in a train-track has an edge parallel to a fixed unit 
vector. We assume that each train-track is extended as far as possible in both direction, 
so that it is a bi-infinite path. 

The idea of the proof is to first understand how the free energy /d is changed when 
the isoradial embedding of the graph G is modified by tilting a family of parallel train- 
tracks. If we can compute the free energy in an extreme situation, that is when all 
the rhombi are fiat, i.e with half-angles equal to or tt/2, then we can compute the 
free energy of the initial graph by integrating the variation of the free energy along the 
deformation from the trivial flat embedding back to the initial graph, by tilting all the 
families of train-tracks successively until all rhombi recover their original shape. Let 
^flat trivial flat isoradial graph, and let Q^'^^ be the corresponding Fisher graph. 

Consider a train-track T of G^. All the rhombi of the train-track have a common 
parallel, with an angle equal to a. Let us compute the derivative of fu as we tilt the 
train-track T as well as all its copies. The function \og det K (z , w) is integrable over 
the unit torus, is a differentiable function of a for all {z,w) S \ {1,1}, and the 
derivative is also uniformly integrable on the torus. The dimer free energy is thus 
differentiable with respect to q, and 

dfn ^ ff d " dz dw 



If — —log det K(z,w) 
JJj2 da 



da 2 ././t2 da 2i7:z 2mw 

dK(z,w)u,v dz dw 



2 ^ j /t2 ' da 2mz 2i'KW 

«,fev(gi)' ■ 

idKu,v sr^ ^, .dlogue 



da ^ da 

e=uv£E{gi) ee-B(Gi) 

In the last line, we used the uniqueness of the inverse Kasteleyn matrix whose coeffi- 
cients decrease at infinity, given by Proposition 5 of [BdTOSj . We also used the fact 
that if e = uv, then = \Ku,v\ and P(e) = Ku,vK~l^. Note that the sum is restricted 
to edges coming from E{Gi), because they are the only edges of Gi with a weight 
depending on an angle. 

There is in this sum a term for every rhombus in the fundamental domain Gi, and 
that term only depends on the half-angle of this rhombus. The variation of the free 
energy of the dimer model along the deformation from Q^'"^^ to Q is thus, up to an 
additive constant, the sum over all rhombi in a fundamental domain of the variation of 
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a function f(^) depending only on the geometry of the rhombus: 



with 



d^ ' ' d9 ' 

where V{6) is the probability of occurrence in a dimer configuration of the edge of Q 
coming from an edge of G, whose half rhombus-angle is 9; and 9^^^ is the angle, equal 
to or ^, of the degenerate rhombus in G^^^ associated to the edge e. 

Recalling that by definition of u, and by Equation ([8]): 



2' ^ ^ 2 27rcos6l' 

and using the fact that ^ log cot | = — , we deduce an explicit formula for the 
derivative of f(^): 

df{9) _ fl _^ TT -29\ 1 _ 1 _^Tr-29 



de \2 2ttcos9 sm.9 2sm9 7rsin26' 



Using integration by parts, a primitive of is given by: 



f(0) = \ log tan ^ + log tan - ^ {l{9) + L (| 



The problem is that f(0) goes to — oo when 9 approaches 0, and thus the free energy 
diverges when the graph G becomes flat. One can instead study sd, the entropy of the 
corresponding dimer model, defined by 

SD = -fD- Yl ^(e)logi^e, (9) 

which behaves well when the embedding degenerates, since it is insensible to gauge 
transformations, and thus provides a better measure of the disorder in the model. 

The sum in ([9]) is over all edges in the fundamental domain of the Fisher graph 
Q. But since edges of the decoration have weight 1, only the edges coming from edges 
of G contribute. Thus, as for Jd, the variation of s^j can be written as the sum of 
contributions s(^e) of all rhombi in a fundamental domain Gi. 

^Sd = Sd{Q)-Sd{Q''^')= Y1 S(0e)-S(0f-*), (10) 

ee£{Gi) 
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where 

s{e) = -f(e)-vi9)\ogi^ie) 

= -2 '"'5""' 2 - — '""^ + ; (^<''> + H 2 - ")) - (2 + ) 2 



2tt \ cos 9 j TT \ V 2 

/ log cot I N / N / \ 

Since lim log cot 6 — = — log 2, and L(0) = L ( ^ ) = 0, we deduce that the 

6i->o V cos / 

values of s{d^^^) are given by the following limits: 

lim s(0) = — - log 2, and lim s(^?) = 0. 

e^o ^ ^ 2 ^ ' e^f 

Let us now evaluate SDiG^^^)- Since the sum of rhombus- angles around a vertex is 27r, 
and since rhombus half-angles are equal to or tt/2, there is in G^^^, around each vertex, 
exactly two rhombi with half-angle 6^'^^ equal to ^ . Let us analyze the dimer model on 
Gn^^: the "long" edges with weight equal to 00 (corresponding to rhombus half-angles 
0) are present in a random dimer configuration with probability 1. The configuration of 
the long edges is thus frozen. The "short" ones with weight equal to 1 (corresponding to 
rhombus half-angles tt/2) are present in a random dimer configuration with probability 

2 ^ TT- 

As noted above, around every vertex of G^^^, there are exactly two "short" edges. This 
implies that "short" edges form a collection of k disjoint cycles covering all vertices of 
the graph G^^^ , for some positive integer k. Such a cycle cannot be trivial, because 
surrounding a face would require an infinite number of flat rhombi, but there is only a 
finite number of them in G^^*. Since these cycles are disjoint, their number is bounded 
by the number of edges crossing the "boundary" of G^^^. Therefore k = 0{n) = o(n^). 

Since all "long" edges are taken in a random dimer configuration, this implies that, for 
every cycle of "short" ones, edges are either all present, or all absent. The logarithm of 
the number of configurations for these cycles grows slower than O(n^), and therefore 
does not contribute to the entropy. The main contribution comes from the decorations, 
which have two configurations each: 

SD(a^"*) = |^i|log2. 

The total number of "short" edges is \ Vn\ (two halves per vertex of Gn), and thus the 
number of long ones is \En\ — \Vn\- From ([TO]) , we deduce that: 

»D = (l-Bil+|Vil)i^ + 



eeE{Gi) 



TT-2eeL , logcot%\ 1/,,,, /vr 

log cot ge- n \ +-[L{Oe)+L{- 



27r \ cos Op I IT \ V 2 
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and using ([9]), we deduce Theorem [TOl 



□ 



The free energy per fundamental domain of the critical Z-invariant Ising model, denoted 
//, is defined by: 

fi = - lim -^\ogZ;[, 

n— ►oo n 

where is the partition function of the critical Z-invariant Ising model on the toroidal 
graph Gn- In |Bax86| . Baxter gives an explicit expression for the free energy of Z- 
invariant Ising models (not only critical), by transforming the graph G with star-triangle 
transformations to make it look like large pieces of glued together, and making use of 
the celebrated computation of Onsager [Qns44j on 1?. As a corollary to Theorem IIOI 
we obtain an alternative proof of this formula at the critical point. 



Theorem 11 ( |Bax86j ). 

log 2 



1 



logtan Be + -(L{9e)+L(- 



IT 



vr 



vr 



Proof. Using the high temperature expansion, and Fisher's correspondence, we have 
from Equation ([2]) that: 



n sinhj(^e))^;:= n sinhj(^e) zi:. 



e€E{Gi) 



As a consequence. 



Moreover, 



fi = fo- J2 sinh J {Be). 

eeE(Gi) 



sinh J (6) = sinh log 



1 + sin ( 

cos 9 



tan 2 tan 9 



so that by Theorem llOt 



\Ei\ + \Vi 



log 2 



+ E 

eeE{Gi) 

, ^ , log 2 
+ l^i| 



7r-29. 



1 



log tan 6*6 - ^ log cot - ^ {L{9e) + ^ {i^ 



2 ^ 2 

eGB(Gi) 



log tan + log tan B^. 



log 2 



E 

eeE(Gi) 



-\ogiaxi9e + -(L{9e) + L(^ 
vr vr \ \ z 



+ 



□ 
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The critical Laplacian matrix on G is defined in [Ken02j by: 

tan 9uv if u V 

- tan Ouu' ifu = v 

otherwise. 



The characteristic polynomial of the critical Laplacian on G, denoted Pj\{z,w), is de- 
fined by P/\{z,'w) = det A(z,ti;), where A{z,'w) is the Fourier transform of the critical 
Laplacian A. 

In a similar way, the characteristic polynomial of the critical dimer model on the Fisher 
graph G, denoted by P{z,w) is defined by P{z,w) = det K{z,w), where K{z,w) is the 
Fourier transform of the Kasteleyn matrix K of the graph G with critical weights. 

In Theorem 8 of |BdT08j . we prove that both characteristic polynomials are equal up 
to a non zero multiplicative constant. Using Theorem 110^ and [Ken02] . we can now 
determine this constant explicitly. 

Corollary 12. The characteristic polynomials of the critical dimer model on the Fisher 
graph Q, and of the critical Laplacian on the graph G, are related by the following explicit 
multiplicative constant. 

yiz,w)eC^ P{z,w) = 2\^'\ H [cot''^-l^PA{z,w). 

Proof. Let c 7^ 0, be the constant of proportionality given by Theorem 8 of [BdTOS] : 

P{z,w) = cPi^{z,w). 

Taking logarithm on both sides, and integrating with respect to (z, w) over the unit 
torus yields: 

\ogc = -2fD- [[ logPAiz,w)^^. (11) 
Moreover, by |Ken02j we have: 



//t2 ^ ' ^2i'Kz2i'Kw 
Plugging this and Theorem [10] into (jlip . we obtain: 



^logtan0e+-(L(0e)+^fJ 
vr vr V V 2 



logC= (l^ll + |yi|)log2+ (log cot y + log cot 



log(2l^il+l^^l W cot I cot 

eeB(Gi) 



iog(2i^^i n H y-i 

ee-B(Gi) 



□ 
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6 Proof of Theorem [5] 



Let us recall the setting: G is an infinite isoradial graph, and Q is the corresponding 
Fisher graph whose edges are assigned the dimer critical weight function; K denotes 
the infinite Kasteleyn matrix of the critical dimer model on the graph Q, and defines 
an operator acting on functions of the vertices of Q. We now prove Theorem [Sj i.e. we 
show that KK~^ = Id, where K~'^ is given by: 

Vx, y € V{g), K-l = MX)fy{-X) Exp,_y(A) log AdA + C,,y. 

This section is organized as follows. In Section 16.11 we prove that the function, 
/a;(A) Expx^y(A), seen as a function of x G ViG), is in the kernel of the Kasteleyn 
operator K. Then, in Section 16.21 we give the general idea of the argument, inspired 
from |Ken02j . used to prove KK"^ = Id. This motivates the delicate part of the 
proof. Indeed, for this argument to run through, it is not enough to have the contour 
of integration Cx^y to be defined as a simple closed curve containing all poles of the 
integrand, and avoiding the half-line dx^y] we need it to avoid an angular sector Sx,y, 
which contains the half-line dx^y and avoids all poles of the integrand. These angular 
sectors are then defined in Section 16.31 The delicate part of the proof, which strongly 
relies on the definition of the angular sectors, is given in Section f6.4[ 



6.1 Kernel of K 



In this section we prove that the function /^.(A) Exp^ y(A), seen as a function of 
X G y{Q)i is in the kernel of the Kasteleyn operator K. In other words. 

Proposition 13. Letx, y he two vertices ofQ, and let xi, X2, x^ be the three neighbors 
of X in Q, then for every A G C.- 



^i^x,xJx,(A)Exp,^,y(A) = 0. 



i=l 

Proof. There are three cases to consider, depending on whether the vertex x is of type 
'w', 'z' or 'u'. 

If X = Wfc(x) for some A; G {1, . . . , (i(x)}. Then the three neighbors of x are xi = ^^(x), 
X2 = -Zfc+ilx) and X3 = Vk{^)- Since x,xi,X2,X3 all belong to the same decoration, we 
omit the argument x. By our choice of Kasteleyn orientation, we have: 

By definition of the angles in M/47rZ associated to vertices of see ([3]), we have: 

Oi = J ^"^^ if ^Wfc.^fc+i ~ ^ 

I a»«fe + 2vr if e«,,,^,+i = -1. 
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Using the definition of the function /, we deduce that: 

K^„z,^,U^,W = -U,W. (12) 

As a consequence, 



3 



1=1 

A, (A)] Exp,,y(A) 

= [-/..(A) - + fv,{\)] Exp,,y(A). (13) 

Since by definition we have fv^W = fwkW + /zfc(A), we deduce that (fT3]) is equal to 0. 

If X = Zk{x) for some k G {1, . . . , d{x)}. Then the three neighbors of a; are xi = tt;fc_i(x), 
X2 = Wki^) and X3 = Vk{x). By our choice of Kasteleyn orientation, we have: 

Similarly to the case where x = Wk, we have: 

-^Zfe,-u)fe_i/»i)fe_i(A) = /^^(A). 

As a consequence, 

3 



X] ^^:^«/^.(^) ExPx„y(A) 



/».-i(A) + /«-.(A) + A, (A)] Exp,,y(A) 

= [/..(A) + /«„(A) - /.,(A)] Exp,,y(A). (14) 

Again, we deduce that is equal to 0. 

If X = Vk{x) for some k £ {!,••• ,d(x)}. Then the three neighbors of x are xi = ^^(x), 
2^2 = Wfc(x) and 2:3 = ?;^(x'), where i and x' are such that '^^(x) ~ Vi{x'). Although 
all vertices do not belong to the same decoration, we omit the arguments x and x', 
knowing that the index k (resp. i) refers to the vertex x (resp. x'). By our choice of 
Kasteleyn orientation, we have: 

K -IK - -1 K -£ cot 



Using the fact that Expx' y(A) = Exp^/ x(A) Exp^ y(A), we deduce: 

3 

^i^x,xJx,(A)Exp,^,y(A) = 

i=l 

= [K^^^^J^^{\) + K^^^yjJ^^{\) + K^^^^^ Exp^,_^(A)^(A)] Exp^_y(A) 



^(A) - /^,(A) + cot [ \ " ) Exp^, ,,(A)^(A) Exp^^y(A) 
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Moreover, by definition of tlie angles in M/47rZ associated to vertices of Q, see (jH), we 
have: 

and using the definition of the function /, we deduce: 

. _ 2 2 \ _ / e' 2 2 \ 

so that: 

-^i)fe,j;<, Expx/^x('^)/f« (^) = 



i cot 



4 / (e*"»fe - A) (e*"-fc - A) \^e^"™fc + A e*"-^ +Xj 

= U,{X)-U{X). (15) 

As a consequence, 
3 

5^i^x.,.,/x,(A)Exp,^,y(A) = [/,,(A) -/^,(A) + ^,(A) -/,,(A)]Exp,,y(A) = 0. (16) 

i=l 

□ 



6.2 General idea of the argument 

The general argument used to prove Theorem [5l i.e. KK~^ = Id, is inspired from 
[Ken02) , where Kenyon computes a local explicit expression for the inverse of the Kaste- 
leyn matrix of the critical dimer model on a bipartite, isoradial graph. It cannot be 
applied as such to our case of the critical dimer model on the Fisher graph Q of G, 
which is not isoradial, but it is nevertheless useful to sketch the main ideas. 

Let X, y be two vertices of Q, and let us assume that the contour of integration Cx^y of the 
integral term of K~^., is defined to be a simple closed curve oriented counterclockwise, 
containing all poles of the integrand and avoiding an angular sector Sx,y, which contains 
the half-line d^^^, see Figure [9] below. 

The argument runs as follows. Let xi,X2,x^ be the three neighbors of x in Q. When 
x / y, the goal is to show that the intersection of the three sectors n^=i ^Xi,y is non 
empty. Then, the three contours Cx^^y can be continuously deformed to a common 
contour C without meeting any pole, and 

dA 



l^i^x,.. i /.,(A)/,(-A)Exp,^,y(A)logA- 

1=1 ''C^i,y ^ ' 



(^j;K.,.J.,(A)Exp,^,y(A)j /,(-A)logA 



^ = 0, (17) 



(2vr) 
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Figure 9: Definition of tlie contour of integration Cx,y of the integral term of K^ y. The 
poles of the integrand are the thick points. 

since by Proposition \T3\ the sum in brackets is zero. Note that in |Ken02j , Kenyon 
has a result similar to Proposition [T3l giving functions that are in the kernel of the 
Kasteleyn matrix of the dimer model on a bipartite, isoradial graph. 

When X = y, the goal is to show that the intersection of the three sectors is empty, and 
explicitly compute: 

Y.^-^-rj /.,(A)/,(-A)Exp,^,,(A)logA-|^ = 1. (18) 

In our case things turn out to be more complicated, since we cannot define angular 
sectors Sx,y so that (fT7|l and (fT8]l hold as such. We define them in such a way that 
when X and y do not belong to the same triangle of a decoration x, the intersection of 
the three sectors ni=i ^Xi,y is non empty, and (fTTl) holds. Then, in order for {KK~'^)x^y 
to be equal to 5x,y in all other cases, we need to have an additional constant Cx,y in the 
expression of the inverse K~l^. Let us mention that angular sectors are not constructed 
explicitly in jKen02j . since geometric considerations suffice. The construction of Case 
1 of Section 16.3.31 could also be used in that setting. 

6.3 Definition of the angular sectors Sx^y 

Let X, y be two vertices of Q. In this section, we construct the angular sector Sx^y 
containing the half-line dx^y and avoiding all poles of the integrand of K~^. In order to 
do this, we first recall some facts about isoradial graphs, then we encode the poles of 
the integrand of in an edge-path ^x,y of the diamond graph C . Using this path 
lx,y, we define the angular sector Sx^y. 

6.3.1 Train-tracks and minimal paths 

Let G be an infinite isoradial graph, and be the associated diamond graph. Let 
us first recall the following definition of |KS05j . A train-track of is a path of edge- 
adjacent rhombi of G*, which does not turn: on entering a face, it exits along the 
opposite edge. As a consequence, each rhombus in a train-track has an edge parallel 
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to a fixed unit vector, called the direction of the train-track. We assume that each 
train-track is extended as far as possible in both direction, so that it is a bi-infinite 
path. By |KS05j . we know that a train-track can not intersect itself, and that two 
distinct train-tracks cross each other at most once. 

Let X, y be two vertices of G, we now define the notion of minimal path of from x 
to y. We say that a train-track of separates x from y, if when deleting it, x and y 
are in two distinct connected components. Now, observe that each edge of G^ belongs 
to a unique train-track whose direction is given by that edge, it is called the train-track 
associated to the edge. An edge-path 7 of G^, from x to y, is called minimal, if all the 
train-tracks associated to edges of 7 separate x from y, and all those train-tracks are 
distinct. In general there is not uniqueness of the minimal path, but all minimal paths 
consist of the same steps, taken in a different order. 

6.3.2 Encoding the poles of the integrand of K~y 

Let x,y be two vertices of Q, and let x, y be the corresponding vertices of G. In this 
section, we define an edge-path ^jx^y of the diamond graph G^ encoding the poles of the 
integrand /^(A)/y(-A) Exp^ y(A) log A of K"^. 

By definition of the exponential function and of a minimal path, the poles of Exp,;. y(A) 
are encoded in the steps of a minimal path of G^, oriented from y to x. It would be 
natural to see the poles of the integrand of K^^y, as representing the steps of a path of 
G^ passing through y and x, obtained by adding the steps corresponding to the poles of 
/a;(A) and fy{—X). However the situation is a bit trickier due to possible cancellations 
of these poles with factors appearing in the numerator of Exp^ y(A). This is taken into 
account in the following way. 

By definition, the function /a;(A) has either 1 or 2 poles. Let us suppose it has 2, 
denoted by e*"^ , e*"^ (the case where it has 1 is similar and easier to handle) , denote 
by T^, the two associated train-tracks, and let = {Tx,Tx}. Similarly, let us 
suppose that fy{—X) has 2 poles — e*"i, — e*"2, defining two train-tracks Ty, Ty, and 
let Ty = {Ty\T^}. 

Let 'yx,y be a minimal path from y to x. If separates y from x, then the pole e*"^ 
is cancelled by the exponential, and we leave 'yx,y unchanged. If not, this pole remains, 
and we extend 'yx,y by adding the step e*"^ so that the initial vertex of 6*°^^ is incident 
to X, let X be the ending vertex of e'"^. This yields a new path of G^, also denoted 
'yx,y, from y to x, passing through y and x. We then do the same procedure for the 
pole e*°2. 

When dealing with a pole of A), one needs to be careful since, even when the 
corresponding train-track separates y from x, the exponential function might not cancel 
the pole if it has already canceled the same pole of /x(A). This happens when Tx and 
Ty have a common train-track. The procedure to extend ''yx,y runs as follows. If Ty 
separates y from x, and TyCiTx = 0, then the pole — e*"'i is canceled by the exponential 
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function, and we leave 'jx^y unchanged. If not this pole remains, and we extend ^x,y by 
adding the step —e^°"^ so that the ending vertex of —e^""^ is y, let y be the initial vertex 
of — e*"i. The same procedure is done for — e*"2, yielding a new path ^x,y of G*, from 
y to X, passing through y and x, and encoding the poles of the integrand of K~^. 

6.3.3 Obtaining an angular sector Sx^y from the path ^x,y 

Let x,y be two vertices of Q. In this section, we construct an angular sector Sx,y, 
containing the half-line dx^y, and avoiding all poles of the integrand of K~^, from the 
path ^x,y encoding these poles. The main step consists in constructing a convex polygon 
Px^y from 7x,j/. The angular sector Sx,y is then naturally defined from Px,y 

In order to construct the convex polygon, let us recall the following definition of |Mer04| 
and lemma of |dT07bj . A finite, simply connected sub-graph of is called train-track- 
convex, if every train-track of which intersects this sub-graph crosses its boundary 
twice exactly. 

Lemma 14 ( |dT07bj ). Any finite, simply connected, train-track convex sub-graph Q of 
can he completed by a finite number of rhombi in order to become a convex polygon 
P whose opposite sides are parallel. 

The construction of the convex polygon P from Q is an explicit algorithm which consists 
in adding a rhombus along the boundary each time the polygon is not convex, see 
Figures [TU] and [T2] for examples. The proof consists in showing that it is possible to 
add a rhombus at every step, and that the algorithm ends in a finite number of steps. 
As a consequence of the construction, we have the following facts: 

1. The train-tracks of associated to boundary edges of Q, and those of P asso- 
ciated to boundary edges of Q have the same parallel direction. 

2. Each train-track corresponding to an edge of P crosses the boundary of P twice 
exactly. 

3. The boundary of the convex polygon P is independent of the order in which the 
rhombi are added. 

We now construct a simply-connected, train-track convex subgraph Qx,y from the path 
7a;, J/, encoding the poles of the integrand of K^y. Applying Lemma E] yields a convex 
polygon Px,y, from which we define the angular sector Sx,y 

Recall that Tx (resp. Ty) consists of the one/two train-track(s) associated to poles of 
/a;(A) (resp. fy{—X)). There are three cases to consider, depending on the intersection 
properties of Tx and Ty . 

Case 1: Tx and Ty have train-track in common. Then, by construction, the path 'jx,y 
of Section 16.3.21 is a minimal path of G^ from y to x, passing through y and x. As such, 
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it is a simply-connected, train-track convex sub-graph of G^. Indeed, ^x,y is a subgraph 
of G*, whose interior is the empty set, and whose boundary is ^x,y itself, consisting of 
"doubled" edges enclosing the empty interior. Since ^x,y is minimal, all train-tracks 
associated to edges of '^x,y are distinct, so that they all cross the path ^x,y once, i.e. 
the boundary of ^x,y twice. We define Qx,y to be '^x,y, and apply Lemma [HI in order 
to obtain a convex polyg on Px,yj see Figure [lO] (left), the numbers indicate the order 
in which the rhombi are added. 



Figure 10: Left: construction of the convex polygon Px^y from a minimal path 7a;, y 
Right: corresponding angular sector Sx,y 

Using 1, 2, 3, and the fact that 'jx^y is minimal, we deduce that the vertices y and x are 
on the boundary of Px,y, and that the boundary of Px^y is independent of the choice 
minimal path of from y to x. Moreover, we deduce that each edge of the boundary 
of Px,y from y clockwise to x is parallel to a unique edge of 7x,j/- The same holds for 
the boundary of Px^y from x clockwise to y. 

Let 

ci) ■ ■ ■ ) Cn be the edges on the boundary of Px,y from y to x, oriented clockwise. 
Then these edges encode the poles of the integrand /x(A)/j;(— A) Exp^ y (A) log A of K~y. 
Since Px^y is convex, they are all included in an angular sector of size smaller than vr. 
We define Sx^y to be an angular sector strictly included in the complement of this sector, 
see Figure [TUl (right) where the sector Sx^y is given in dashed lines. 

Case 2: Tx and Ty have 1 train-track T in common. We specify the construction of the 
path 'jx^y as follows. Let us first consider the poles of fx{^) and fy{—X) which are not in 
the common train-track T (if any), then there is a minimal path from a vertex yo to a 
vertex xq, passing through y and x, encoding these poles and those of the exponential 
function. We then consider two cases: the common train-track T may or may not 
separate y from x. 

Suppose that the common train-track T does not separate y from x, see also Figure 
[TTl Then the vertices yo, y, xq, x are on the same side of T, and the minimal path 
from yo to xq follows the boundary of the train-track T. Moreover the pole of /^(A) 
corresponding to T is e*", that of fy{—X) is the opposite — e*", and the exponential 
function cancels none of these poles. As a consequence the path 'jx^y is obtained by 




Y 

■ r 



33 



adding the vector — e*" at the beginning of the path, and the vector at the end, see 
Figure [TT] (thick hnes). 



Figure 11: The common train-track T does not separate y from x. The path ^x,y is 
given in thick hne. 

Suppose that the common train-track T separates y from x, see also Figure [T^l Then 
the vertices yo, y are on one side of the train-track T, and the vertices xq, x are on the 
other. The minimal path from yo to xo follows the boundary of the train-track T for 
a number of steps (which can be 0), then crosses it, and follows the boundary on the 
other side until xo. Moreover, /a;(A) and fy{—X) have the same pole e*" corresponding 
to T. The exponential function cancels one of the two poles, say the one of /x(A), 
so that the pole of fy{—X) remains. As a consequence, the path jx,y is obtained by 
adding the pole e*" at the beginning of the minimal path from yo to xq, and leaving 
Xq unchanged. The path 'yx,y thus crosses the train-track T from y to yo, follows the 
boundary for a number of steps, then crosses the train-track T again, and follows the 
boundary until x = xo. Just by changing the order of the steps, we can choose 'yx,y as 
in Figure [12] below (thick lines). 



Figure 12: The common train-track T separates y from x. The path jx,y is given in 
thick line. 




y 



y 
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Note that in both cases, the path ^x,y starts with a parallel direction of the common 
train-track T from y to yo, then follows its boundary, and ends with the opposite 
parallel direction until x. Define Qx,y to be the subgraph of made of the rhombi 
of the train-track T that are bounded by ^x,y Then, Qx,y is train-track convex and 
we can apply Lemma [T3] in order to obtain a convex polygon Px,y, see Figure [13] below 
(left), the numbers indicate the order in which the rhombi are added. 




Figure 13: Left: construction of the convex polygon Px^y from the train-track-convex 
subgraph Qx,y Right: corresponding angular sector Sx,y 

Using 1, 2, 3, the fact that the part of ^x,y from yo to x is minimal, and the fact that the 
path '^x,y crosses the train-track T twice, we deduce that the vertices y, x, yo, xo are on 
the boundary of Px,y Moreover, we deduce that each edge on the boundary of Px^y from 
y to X passing through yo is parallel to a unique edge of '^x,y Let ei, • • • , denote this 
set of edges. Then, they encode the poles of the integrand /x(A)/y(— A) Expx y(A) log A 
of K~^, and define a sector of size vr. We define Sx,y to be an angular sector strictly 
included in the complement of this sector, see Figure [13] (right) where the sector is 
given in dashed line. 

Note that the above definition of angular sector Sx^y is ambiguous when the path ^x,y 
crosses the train-track T twice, without following its boundary in between. This can 
only occur in the following four cases. Since all vertices involved belong to the same 
decoration, we omit the argument x in the notations. Recall that e*°"fc = e*°^*=+i . 

{wki Zk+i) then, ^x,y is the path — e*"^fc+i , e*°™fe , 
{zk+i,Wk) then '-^x,y is the path -e*""''^ , e'^^^^+i , 

[x,y) = < 

{wk,Wk) then ^x,y is the path — e*"""* , e*"™*: , 
{zk+i, Zk+i) then 'yx,y is the path — e*"^'=+i , e*"^*+i . 

In these four cases, we set the conventions given by Figure [E] below for the angular 
sector Sx,y These are natural in the light of the proof of Theorem [^ see Section [6.4i 
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Figure 14: Definition of the angular sector Sx,y (daslied line) in the four ambiguous in- 
stances of Case 2: from left to right, = (wfc, Zfc+i), (zfc+i, Wfc), (wfc, lUfc), Zfc+i). 

Case 3: and Ty have 2 train-tracks in common. This can only occur in two cases: 
either the vertices x and y are neighbors in Q and the vertices x, y are neighbors in 
G, meaning that x = ffc(x) and y = v^(y), with k and i such that ffc(x) ~ ve{'y); or 
X = y = Vfc(x). 

Suppose that {x,y) = (u^ (x) , (y ) ) . Then /^j^(x)(A) and /„^(y)(-A) have the same 
poles e*"'"*:^''', e*"^fc''''. The exponential function Expxy(A), cancels one of the pairs 
of poles, and adds two new ones —e'^^k^^)^ _g*°zfe{x)_ Thus the path ^x,y follows the 
boundary of the rhombus associated to the edge xy (note that this rhombus is already 
our convex polygon Px,y)- By convention, we choose the order of the edges of the path 
'yx,y as in Figure [T5] (left). This is natural in the light of the proof of Theorem [51 see 
Section 16.41 This defines a non-empty angular sector Sx,y, containing no pole of the 
integrand of K~y, see Figure \TE\ (right). 




Figure 15: Left: Choice of path 'Jvk{x),ve{y)- Right: corresponding angular sec- 
tor S^,.(x),^,(y)- 

Suppose that {x,y) = {vk{x),Vk{x)). Then /t,^(x)(A) has poles e*"™fcW, e*""'fc<''', and 
/t,j.(x)(— A) has opposite poles. Moreover, Expxy(A) = 1, so that it cancels no pole. 
Thus the path 'jx^y follows the boundary of the rhombus defined by the vectors e*"™*:'"' 
and e°'''k(^) . By convention, we choose the order of the edges of the path "jx^y as in 
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Figure [TBI (left). This defines a non-empty angular sector Sx^y, containing no poles of 
the integrand of K~^, see Figure [T6] (right). 




x=y 




'k,Vk 



Figure 16: Left: Choice of path "yvf,{x),Vi.{y)- Right: corresponding angular sec- 
tor S^,.(x),„,(x)- 

Remark 15. By the above definition, in all cases but the four exceptions to Case 2 and 
Case 3, the angular sector Sx,y always contains the half-line dx,y starting from 0, in the 
direction from x to y. Thus, we can define the contour of integration Cx,y to be a closed 
contour oriented counterclockwise and avoiding the closed half- line dx^y 



6.4 Proof of Theorem [5] 

In this section, we prove Theorem [5l i.e. KK^^ = Id. Let / and C be the infinite 
matrices whose coefficients are the integral and constant part of K^^ respectively: 

Ix,y = jf fx{\)fy{->^) Exp^ y(A) log AdA, 

^x,y — Cx^y 

Our goal is to show that {KK~^)x,y = {KI)x,y + {KC)x,y = Sx,y Let xi, X2, be the 
three neighbors of x in ^. Then, 



iKI)x,y = Y,i i^.,.J,,(A)/,(-A)Exp,^,y(A)logA 



dA 



2 ■ 



{KC)x,y = Y.Kx,xfix,,y (19) 



i=l 



In Proposition 1171 we handle the part {KI)x,y'- we prove that as soon as x and y do 
not belong to the same triangle of a decoration, n?=i ^Xi,y / 0) so that by the general 
argument of Section 16.21 (Kl)x^y = 0; when x and y belong to the same triangle of 
a decoration, we explicitly compute {KI)x,y Then in Lemma [TSl we handle the part 
{KC)x,y. we show that the constants Cx,y are defined so that {KK~^)x,y = Sx,y The 
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proof of these results is always complicated by the fact that x can be of three types, 
'it;','-?' or i.e x = u)fc(x), Zk{x) or Wfc(x) for some k G {1, • ' ' yd{'x)}. The next 
proposition gives relations between these three cases for KI. Note that whenever no 
confusion occurs, we omit the argument x. 

Proposition 16. For every vertex y of Q, the quantities {KI)wf.,y, {KI)z^^y and 
{KI)y^^y satisfy the following: 



1. {KI) 



Wk,y 



+ 



/..(A) + (- 



+ 



fwkW /j,(-A)Exp (A)logA 



2. (KI) 



■"k,y 



fzkW + 



+ 



fwkW /y(-A)Exp^y(A)logA 



dA 
dA 



(27r)2 



3. // Su.fc,y n Sy^^y / 0, then {KI) 



Wk,y 



-{KI) 



Zk,y 



<KI).„y 



+ 



fzkW /y(-A)Exp^y(A)logA-^^. 



4. // s^^^y n s^^^y i= 0, then {KI) 



Wk,y 



iKl),^,y={Kl) 



Vk,y 



+ 



fwkW /y(-A)Exp^y(A)logA 



dA 



5. //s^fc.s/ n s^,^,y n / 0, then {KI)^^^y = {KI)-,^^y = {KI)^^^y 



0. 



Proof. Points 3, 4, 5 are a consequence of 1 and 2, and the general argument of Section 
I6.2i Let us prove 1 and 2. In equations p^ . (fH|) . (fT5|) . ([16]) of the proof of Proposition 
[T3| we explicitly computed 



j=i 



for X = Wk{x), -za:(x) and Vk{x), respectively. Using this. Equation p^ . and the fact 
that fvkW = fwkW + fzkW, we obtain the following. 

If X = w/c(x), then the three neighbors of x are xi = -Zfc(x), X2 = Zk+i{'x.), X3 = ffc(x), 



(20) 
dA 



and by ()T3|1 we have: 

{KI)wk,y = 



(- i fzkW - i UkW + f [/^„,(A) + /,,(A)])/,(-A)Exp,,y(A)logA 



(27r)2 
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If X = Zk{x), then the three neighbors of x are xi = Wk-i{x),X2 = Wk{x),X3 = fA;(x), 
and by (fTll) we have: 



{KI).,,y = (21) 

{ /..(A) + f U,{X) - i [/^,(A) + /,,(A)])/,(-A)Exp,,y(A)logA-|^. 

If X = Vk{x), then the three neighbors of x are xi = Zfc(x), 2:2 = Wk{x),xs = Vi{x'), 
where i and x' are such that fA;(x) ~ ve{x.') in ^. Recalhng that the index £ refers to 
the decoration x', we also omit the arguments x and x'. Using (llSp and (I16p . we have: 

iKlU,y = (22) 

dA 



U fz,{\) - f f^,W + f [/^,(A) - /,,(A)])/,(-A)Exp,,y(A)logA 



(2^) 



2 • 



As a consequence of (pOj) . (f2T]) . (p2|) . 1 and 2 of Proposition [16] are proved, if we show 
that, for every vertex y of ^, 

which is equivalent to proving that Sz^j^^^y = Swk,y^ aiid s^^.^^ = s^g^y Recall that the 
construction of the sector Sx,y, given in Section 16.31 relies on the path 'yx,y encoding 
the poles of the integrand of Kx,y- Recall also that, in all cases but the four exceptions 
to Case 2, and Case 3 of Section [6.3.31 when two pairs of vertices of Q have the same 
path 'yx,y, then they have the same sector Sx^y 

Proof of Sz^^^^y = Swk,y The vertices z^+i and Wk belong to the the same decoration 
X, so that the exponential function has the same poles in both cases. Moreover, the 
functions fzk+ii'^) ™d fw^W have the same pole e^^^'^+i = e*"™* . As a consequence, 
by construction, the sectors Sz,.^-^,y and Swf,,y are the same in all cases but the four 
exceptions to Case 2, which have to be checked 'manually'. Note that Case 3 cannot 
occur when x is of type 'zu' or 'z'. Referring to Figure [El we see that by definition: 

Proof of Syj,^y = Sy^^y. Thc vertex (resp. v^) belongs to the decoration x (resp. x'), 
and k,i,x,x.', are such that Ufc(x) ~ Vi{x') in Q. Let us prove that the paths 7^^(x),j/ 
and 7t,^(x').j/5 encoding the poles of the integrand, are the same. The function /t,^(x)(A) 
has poles e'^^k(^\ e°-''k(.^\ and the function /„^(x')(A) has opposite poles e*"™^''''' = 
_g*"ii.fe(x)^ g*"z^(x') _ _g*"zfc(x)^ consequence the vertices Vfc(x) and f^(x') define 

the same two train-tracks denoted by Tx = {r^,T^}. Let Ty be the one/two train- 
track(s) associated to the vertex y. Since the construction of the path ^x,y is split 
according to the intersection properties of Tx and Ty, we use the same decomposition 
here. 
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IfTx and Ty have train-track in common. Then the paths 7t,j.(x),y and 'yvi{x'),y 
minimal, and are constructed according to Case 1 of Section r6.3.3l Recall that in Case 1, 
the construction of the angular sector Sx,y from the minimal path 'j^^y is independent 
of the choice of minimal path of from y (the initial vertex) to x (the final vertex). 
Moreover, by construction, the initial vertex of 7^j.(x),j/ and of 7^^(x'),j/ is the same, let 
us denote it by y. So we are left with proving that the ending vertex of both paths is 
the same. 

• IfT^, do not separate y from x. Refer to Figure [171 Then, Exp^ y(A) does not 
cancel the poles of fvf,{x)Wi and the path 7t,j,(x),y is obtained by concatenating 
a minimal path from y to x, and the vectors 6*"'"'='''', 6*°^*=^''^ The ending vertex 
is Vfc(x) = x + e*""'fe('') + e^^^ki^) = x'. When the train-tracks and T| do not 
separate y from x, then they separate y from x'. As a consequence, Exp^/ y(A) 
cancels both poles of /^^(x')(A), and the path 7t,^(x'),y is a minimal path from y 
to x'. 




Figure 17: and Ty have train-track in common. Construction of the paths 7t,^(x),y 
and 7t,^(x'),j/) when T^, do not separate y from x. 

• // only separates y from x. Refer to Figure [T8l Then, Exp^ y(A) cancels 
the pole e*"'"fc(''' of /„^(x)(A), and the pole e*"^^^''' remains. The path 7^j.(x),j/ is 
obtained by concatenating a minimal path from y to x and the vector e^'^^k '''^ . 
The ending vertex is i'fc(x) = x -|- e*"'^*^'''. When only separates y from x, 
then only separates y from x'. As a consequence, Expx/_y(A) cancels the 
pole — e*"^fcW of /t,^(x')(A), and the pole — e*""'*'''' remains. The path 7^^(x'),y is 
obtained by concatenating a minimal path from y to x', and the vector — e^^'^fe^"' . 
The ending vertex is, ^^(x') = x' — e*"'"*'''' = x -|- e*"^fe^''^ = Vk{x). 




X 



Figure 18: and Ty have train-track in common. Construction of the paths 7„^(x), 
and 7^^(x'),y) when only T^ separates y from x. 
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All other cases are symmetric. 

IfTx and Ty have 1 train-track T in common. Then 7^^.(x),j/ ^'^d 7t,^(x'),y paths con- 
structed according to Case 2 of Section 16.3.31 Recall that in Case 2, the path jx,y starts 
from the initial vertex y with a parallel direction of the common train-track T, then 
follows its boundary for a positive number of steps, and ends with the opposite direc- 
tion up to the ending vertex x. By construction, the initial vertex of 7t,j.(x),3/ and of 
lvi{-x.'),y is the same, let us denote it by y. We are left with showing that both paths 
have the same ending vertex. The argument is similar to the previous case. We do not 
repeat it here, but only provide Figures [T9] and [20] below, which illustrate it. 




Figure 19: and Ty have 1 train-track in common. Construction of the paths 7t,^(x).. 
and 7t,^(x'),j/) when T^, do not separate y from x. 




Figure 20: and Ty have 1 train-track in common. Construction of the paths 7^^(x),j/ 
and 7t,^(x'),j/j when only T^ separates y from x. 

Note that the specific instances of Case 2 do not occur, since x is of type 

IfTx and Ty have 2 train-tracks in common. This can only occur when y = t'fc(x) or 
y = Vi(x'). So we need to check that: 

Referring to Figures [15] and [T6] we see that this is true by definition. □ 

We now state Proposition [T7] computing the matrix product KI of Equation p9|) . 
Whenever no confusion occurs, we drop the argument x in Wfc(x), Zk{:x.), ■Vfc(x). 
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Proposition 17. For all vertices x and y of Q, we have: 

if {x,y) = {wk,Wk), {zk,Zk), ivk,Wk), {vk,Zk) 
if {x,y) = {zk,Wk), {■Wk,Zk) 
if {x,y) = {vk,Vk) 
else. 

Proof. Let X be a vertex of ^. Then x = ■Wfc (x) , (x) or (x) for some A; G {I,-"" ,d(x)}. 
Let us show that whenever y ^ {wk{x), ^^(x), ^^(x)}, then the intersection of the three 
angular sectors Sy^^(^)^y n S2^(x),j/ H s^j,(x),3/ is non empty. Using 5 Proposition [HI this 
implies that {KI)x,y = 0. 

Recall that the construction of the angular sector Sx^y, given in Section [6.3.31 relies on 
the path ^x,y encoding the poles of the integrand of K~^. The three vertices ^^^(x), 
2;fc(x), ffc(x) all belong to the same decoration x, so that the exponential function 
is the same in all three cases. The function /^j.(x)(A) has pole e^°'^k(^\ the function 
fz^.{x)W has pole e*"^fe'''), and the function /t,j.(x)(A) has poles e^^'^fe*"), e^'^^'^i^) , As a 
consequence the vertices Wk{^), Zk{x), Vfc(x) define the same two train-tracks denoted 
by Tx = {T^,T^}. Let Ty be the one/two train-track(s) associated to the vertex y. 
Since the construction of the path 'yx,y is split according to the intersection properties 
of Tx and Ty, we use the same decomposition here. 

IfTx and Ty have train-track in common. Then the paths 7«n.,j/, "yzk,y^ 7«fc,y min- 
imal, and are constructed according to Case 1 of Section 16.3.31 Recall that in Case 
1, the construction of the angular sector Sx,y from the minimal path is independent of 
the choice of minimal path from y (the initial vertex) to x (the ending vertex). More- 
over, returning to the construction of the convex polygon P from a minimal path 7, we 
have that if 71 and 72 are two minimal paths such that 71 C 72 (meaning that edges 
of 71 are also edges of 72), then the boundary of the convex polygon Pi is included 
in the boundary of the convex polygon P2. As a consequence, the angular sector S2 
corresponding to P2 is included in the angular sector si corresponding to Pi. 

• IfT^, T^ do not separate y from x. Then the minimal paths 'ywk,y, Iz^^y^ 7ffe,y 
can be chosen as in Figure [2T] below. 




Figure 21: Tx and Ty have train-track in common. Construction of the paths 7, 
7zfe,2/) lvk,y^ when T^^ and T"^ do not separate y from x. 



{KI) 



x,y 
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Observe that ^Zk,y C 7^^.,^, so that by the above remark s^^^y C Szy,,y Since the 
construction of the angular sector s^^^y is independent of the choice of minimal 
from y to -Ufc, we can also choose ^Vk,y to be the concatenation of the minimal path 
from y to x, and the vectors e*""'fc , e^^-'^ . Then, ^Wf,,y C ^Vk,y> ^-^^d s^^^y C s^j.,?/- 
We deduce that: 

since s^^^^y is an angular sector of size at least vr. 

• // only separates y from x. Then the minimal paths ^w,,,y, "yz^,y, lvk,y can be 
chosen as in Figure [22] below. 




= Wt 



Figure 22: and Ty have train-track in common. Construction of the paths "fw^^y, 
1zk,y, 7vk,y, when only separates y from x. 

We have 'ywk,y C Jzk,y = Iv^^y, SO that s^^^y n Sz^,y Pi s^^^y = s^^^yi^^, since s^^^y 
is an angular sector of size at least vr. 

• If and separate y from x. Then the minimal paths ^w^,y, lzu,y> lvk,y can 
be chosen as in Figure [23l below. 




Figure 23: and Ty have train-track in common. Construction of the paths Jwf^y^ 
lz^,y, 7v^,y, when and separate y from x. 

Then, the three paths are equal, so that s^^^y n Sz^,y H 5^,^.^^ = Sw^.,y 7^ 0- 

IfTx and Ty have 1 train-track T in common. Let us assume that T = T^ = Ty. 
Then, the paths Jzk,y and 'yvk,y are constructed according to Case 2 of Section [6.3.31 
and 7wj.,y is constructed according to Case 1. 

Suppose that y ^ Wk, z^j^i, Zk,Wk-i, so that we are not in the exceptions to Case 2. 
We only handle the case where and do not separate y from x, see Figure [Ml 
the other cases are handled in a similar way. Then, the paths Jw,,,y, lzu,y> Iv^^y are 
constructed as in Figure [Ml below. 
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Figure 24: Tx and Ty have 1 train-track in common. Construction of the paths "fwk,y, 
lz,:,y, 7v,,,y, when and do not separate y from x. 

By construction of the convex polygon associated to these paths, we have: s^^^y = Sz^^y, 
and Sv^^y C Swfe,y We deduce that: 

^'"k,y ^ ^Zk,y ^ ^Wk,y ~ ^■"k,y ^i 
since Sy^^y is an angular sector of size vr — e, for some e > 0. 

Suppose now that y = Wk, Zk+i, Zk, or Wk-i- In these four cases, the angular sectors 
Szk,y, Su,^,y, s^^^y are drawn on Figure [25l 



a) 



b) 



c) 



d) 




Figure 25: Angular sectors Sz^,y, Sw^,y, s^^^^y in the four exceptions to Case 2. 
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From b) and d) we see that when y = z^+i or y = Wk-i, the intersection Sz,.,y n s^^.^^ n 
•^ffc.y 7^ -1^^ t'^^ remaining two cases we do exphcit computations. 

• Computations for y = Wk- In this case, see Fig 

ure[25]a.), we have s^j.^uifel^ '5vj.,uife 7^ ^ 

Hence, by 4 of Proposition I16[ we know that: 



(KI) 



(KI),,,^, = {KI) 



+ 



■^k'^k "^'"k'^k 



fwkW /^,(-A)logA 



dA 



Using the definition of fw^W^ denoting by C a generic simple closed curve 
oriented counterclockwise, containing all poles of the integrand, and avoiding a 
half-line starting from 0, yields: 



1 



(2vr)2 



/«;,(A)/^,(-A)logAdA 



1 



(27r)2 Jc (A - e*"-fc)(^ + e'"™0 



log AdA 



4tt 



[logc(e*"»0-logc(-e*"»0]- (23) 



As a consequence, see Figure [25] a) 
\ogr {e'°"^k ) + log^ 



i ' 

i r 



(_ei"»fe ) + log (e*°»fc ) - logc (-e*"»fc ) 



Computations for y = z^. In this case, see Figure l^c). we have Syjj^^zj^'^Svi^,Zf^ 7^ 
Hence, by 3 of Proposition 1161 we know that: 



{KI) 



Wk,Zk 



-{KI) 



Zk,Zk 



-{KI) 



Vk,Zk 



+ 



fz,{\) /.,(-A)logA 



dA 

W) 



2 ■ 



Using the definition of /^^^(A), yields fzk{^)fzk{-^) 
using Figure [53 c), we obtain: 



e k 



(A-e*"^fe)(A+e'""fe)' 



SO that 



{KI)n.,,z, 

i ' 
" ~4^^" 



= -{KI)z„z, = -{KI),^,z,= 
logc,_ + logc, ^ 



ioz,^ + i(a2j^ - vr) + ioz,^ - i{az^ + vr)] 



(e-^0-logc„^,.J-e^"^0 
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IfTx and Ty have 2 train-tracks in common. This can only occur when y = Vk{x) or 
ve{x'), such that t^fc(x) ~ vi{x'). In these two cases, the angular sectors Szj^^y, Swk,y, 
Sy^^y are drawn on Figure! 



a) 



b) 































/piil|x 















Figure 26: Angular sectors Szj.,^, Sw^,y, Syf,^y in Case 3. 

From b) we see that {KI)u,^,vi = {KI)^^^^^^ = {KI)^^^^^^ = 0. For y = Ufc(x), we do 
explicit computations. 



• Computations for y = v^- We use 1 and 2 of Proposition [TOl Let us first compute: 

+ f )/.,(A)/.,(-A)logA 



(2^) 



2 ■ 



Using the definition of /^^.(A), /^^.(A), and denoting by C a generic contour, we 
have: 



(^jf/..(A)/.,(-A)logAdA 



4^^ 7c \ (A - e*'"-ft ) (A + e^"-fc ) (A - e'^^k ) ( A + e^°-fc ) 

log^(e-^.)_log^(_e--.-) 



47r 



cos 



log A dA 
(24) 
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As a consequence: 



+ 



/..(A)/.,(-A)logA 



dA 

(2^ 



1 

1 

~2' 



-vr H ^7 — — vr 



cos 



cos 



Let us now compute, 



+ 

Exchanging with in (j24p yields: 

/ /«;fc(A)A,(-A)logAdA 



/«'fc(A)/»;fc(-A)logA 



dA 



(2vr) 



2 ■ 



C 



4tt 



logc(e'"-0-logc(-e^°-0 



log^(e-».)_log^(_e-^.) 



cos 



(25) 



As a consequence: 



+ ; /..(A)/..(-A)logA^ 



vr 



cos 



1 

47r 

_ 1 
~ 2' 

Thus, by 1 and 2 of Proposition [HI we deduce: 



+ 7r + 



cos 



□ 

By Proposition 1171 proving that KK~^ = Id, amounts to proving the following lemma 
for KC. 
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Lemma 18. For all vertices x and y of Q, we have: 



i ifix,y) = {wk,Wk), izk,Zk), iwk,Zk), {zk,Wk) 
-\ if {x,y) = {vk,Wk), {vk,Zk) 



else. 

Proof. By definition of C, we have Cx,y = 0, as soon as x or y is of type ^v' . As a 
consequence, using Equation and our choice of Kasteleyn orientation, we deduce: 

{KC)zf,,y = £zk,Wk-iCwk-i,y + Cwk,y 

{KC)vf,,y = Czf,,y — Cwf,,y. (26) 

Let us first prove that for every vertex y of Q, we have £wk,Zk+iCz^_^_i,y = Cwf,,y, or 
equivalently: 

^Zk+i,y ~ ^w,^,Zk+iCwk,y (27) 

When y belongs to a different decoration than x, then both sides of (I27p are equal to 
0. Let us thus suppose that y is in the same decoration as x. If y 7^ -Zfc+i) then by 
definition: 

(_l)n{zk+i,y) I ^ t ^ 

'~^Zk+i,y — ^ ~ 4^ \ ~ ^Wk,Zk+i'^Wk,y! 

SO that (f27|) holds. If y = -Zfc+i, then by definition, C^j._,_i,^j._,_i = —\- Moreover, 
SO that (P7|) also holds in this case. As a consequence. Equation ([^B]) becomes: 

{KC)wf„y = {KC)z^,y = - {KC)v^^y = -Cz^,y + C^fe,!/- 

Let us now end the proof of Lemma [THJ If y does not belong to the same decoration as 
X, or if y belongs to the same decoration as x and is of type 'u', then Cwf,,y = Czf,,y = 0, 
so that: 

{KCU,y = {KC)z„y = {KC),„y = 0. 

If y belongs to the same decoration as x, but not to the same triangle of the decoration, 
then 



(_l)n(wt:,y) 

n — _ ±('_iAn(«'fe,2fe)('_i\n(zfe,y) _ r< —n 

'^Wk,y — ^ ~ 4^ ^ ^ — ^Zk,Wk'^Zk,y — 

since by our choice of Kasteleyn orientation, Szk,Wk = 1- Thus, 

{KCU„y = {KC)z„y = {KC),„y = 0. 



48 



If y = Wk, then by definition (7^^,^.^ = j, and Cz,^,Wk = -Cwk,Zk = thus: 
Finally, if y = z^, then by definition C^„^^,^^^ = \, and C^^^^^^^ = - j, so that: 

□ 

A Probability of occurrence of single edges 

Let us compute the probability of occurrence of single edges in dimer configurations of 
the Fisher graph Q chosen with respect to the Gibbs measure V of Theorem [9l Every 
edge of G is of the form w^Zk, w^Zk+i, w^Vk or VkVe, as represented on Figure [271 below. 
The vertices Zk-, Wk, Zk+i and Vk belong to the decoration of x, and vi belongs to that 
of y. The edge of G joining x and y is the diagonal of a rhombus with half- angle 9. In 
order to simplify notations, let us write 

azf. = a, = /? = a + 29. 




Figure 27: A piece of the Fisher graph Q near a rhombus with half-angle 9, adjacent 
to the decorations of x and y. The sides of the rhombus are represented by the unit 
vectors and e*'^, with [5 — a = 29. 



By Theorem [9l we know that the probability of an edge e = wi' of ^ is given by 



V{e) = Ku,v Pf 



u,v v,v. 



where the coefficient K^^ J of the inverse Kasteleyn matrix is given by Theorem El 
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Probability of the edge w^Zk 



The vertices Wk and Zk belong to the same decoration x, so that there is no contribution 
from the exponential function, since Exp^ x = 1- By our choice of Kasteleyn orientation, 
we have Kz,.^Wk = 1- Using Formula @ and the definition of the function /, yields: 

ViwkZk) = K,„^,K-l,^ = 4^ i U,Wfz,{-X) log(A)dA + ^ 



e 2 



4vr2 L e*^ - A + A 



log(A)dA + ^, 



where Ci = Cw,.,Zk is the contour defined according to Case 1 of Section I6.3.3t repre- 
sented on Figure [28l The integral is evaluated by Cauchy's theorem: 



1 f e'^ log(A) 
i^Jci (A-e*/3)(A + e*") 



dA 



i / logc^(e ^/^)-logg^(-e^«) 



27r 



Therefore, V{wkZk) = \ + ^ — ^ 

4 47r cos d 



-J— {il3 - i{a + vr)) 

4iT cos(^^) 

IT -26 



47r cos 6* 





Figure 28: Left: the contour Ci = Cwf,^Zk involved in the integral term of K^^^^^. Right: 
the contour C2 = Cz^.^;i^,Wk involved in the integral term of K~^^^_^^^^. In both cases, the 
arrows represent the poles of the integrand, and the shaded zone is the angular sector 
avoided by the contour. 



Probability of the edge w^z^+i 

Since the vertices and z^+i belong to the same decoration, there is no contri- 
bution from the exponential function. Moreover, by Equation p2p . we know that 
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Kw^,zf,+Jz^+,W = -/«,fe(A), and by Equation dlZl), we have K^,.,^,^,C^^^^,^^ = C^^.i^, = |. 
Therefore, 



C2 



4^2 



where C2 = Cz^j^^^wi^ is the contour defined according to Case 2 of Section [6.3.31 repre- 
sented on Figure!^ Using Equation (j23|) where we computed this integral for a generic 
contour C, we obtain: 

-4^ jf ^.(A)/«,,(-A) log AdA = ^[\ogc,{e'^) - logc,(-e*^)] = ^{if3 - i{f3 + vr)) = \. 
We conclude that V{wkZk-\-i) = -■ 



Probability of the edges WkV^ an z^Vk 

Again, there is no contribution from the exponential function. By our choice of Kaste- 
leyn orientation, we have K^^^^^k = ~1) by definition we have Cy^^^f, = 0. Using 
Formula ([5j), this yields: 

V{wkVk) = K^,,n„K-l^^ = jf /^,(A)^,(-A)log(A)dA, 

where C3 = Cw,, ,vu is the contour defined according to Case 2 of Section 16.3.31 repre- 
sented on Figure [29l Using Equation ()25p where we computed this integral for a generic 
contour C, we obtain: 



ViwkVk) = — I logc3(e^^) -logc3(-e^^^ 



logC3(e''') -logC3(-e' 



cos , — 2- 



47r "^3' ' ' //3-a 



i A„ ..^ s i/3-i(a + 7r) 
— z/3-z(/3 + ^) 



47r \ cos 6* 

1 71-26 
4 47r cos 6* ' 



By symmetry, this is also the probability of occurrence of the edge z^Vk- 
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Figure 29: Left: the contour C3 = C^^^Vf, involved in the integral term of K^^^^^. Right: 
the contour C4 = Ct,^,!^. involved in the integral term of 



Probability of the edge VkV£ 

By definition, we have Cvf^v^ — 0) by Equation (|15p . we know that: 

Using Formula this yields: 

VivkVi) = K^^^^^K-]^^ 
1 



47r2 
1 

4^ 



Kvk,vJvtWfvk{-^) Expy^x(-^) log AdA 

C4 

[/«,,(A)-/.,(A)]/,,(-A)log(A)dA 
C4 



= 7^/ /«„(A)/.,(-A)log(A)dA--^ / /,,(A)/,,(-A)log(A)dA. 

where C4 = C^^,^^. is the contour defined according to Case 3 of Section [6.3.3l represented 
on Figure [29l Using Equation (f25]l . the first term is equal to: 



^ jf /^,(A)/„,(-A)log(A)dA 



^_^./3) _ logc,(e^/^)-logc,(-e-) 



i ( .r. -la ^ «/3-i(a + vr) 
ip — i\p — vr) 



cos , 2 



47r \ COS0 

1 TT - 261 

+ 



4 47r COS ( 
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By Equation (p^ and Figure [291 the second term is equal to: 
-T^ / /.,(A)/,,(-A)log(A)dA = 

i (, I la., logc,(e'°)-l°gc .(-e"") 
= j;^lo.c.(. )-.o.e.(- ) ^ 

i / — — vr) 

= — I la - i(a + TT) 

47r V cos 6* 

1 vr - 26* 

+ 



4 An cos 6* 

N 1 7r-20 

As a consequence VivkVe) = — I -. 

2 zir cos y 

Let us make a few simple comments about the values of these probabilities. 

1. The value | for the probability of the edge w^Zk+i can be explained as follows. By 
Fisher's correspondence, once the configuration of the edges coming from edges 
of G attached to the decoration of x is fixed, there are two possibilities for the 
dimer covering inside the decoration, which both have the same weight. There 
is always one of the two possibilities containing the edge WkZk+i- Therefore, this 
edge appears in a random dimer configuration half of the time. 

2. Notice that V{wkZk) = ^^{vkVi). This is explained by the fact that the edge w^Zk 
appears only if VkV^ is not present in the dimer configuration, and it appears only 
in one of the two allowed configurations, once the state of the edges coming from 
edges of G is fixed. 

3. Using the two previous points, and the fact that the probability of the edges 
incident to a given vertex must sum to 1, one can deduce the probability of all 
the edges from the probability of the edge VkVf,. 
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